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PREFACE 

The work on the problems in this report was started by the 
author at Douglas Aircraft, Santa Monica, in January, 1946. 

The main bulk of the work was done at the Jet Propulsion Labora- 
tory, California Institute of Technology. 

Many of the methods and results of this report have been 
presented at seminars, in particular at Johns Hopkins University 
and California Institute of Technology. Some of the formulas 
obtained have been used in References 1, 2, and 3. 

The main concepts in the theory of conical wings are drawn 
from the work of A. Busemann (Cf. Ref. 4). Through lectures 
and discussions the author also received many valuable ideas 
from W. D. Hayes, R. T. Jones, and in particular H. J. Stewart. 

Some of the results in Sections I, II, III, and IV were obtained 
independently by the author and others, in particular W. D. 

Hayes, and also by many research workers here and abroad 
working with widely different methods. Sections V and VI are 
believed to be essentially new, as regards both the basic 
solutions and the proposed application of these solutions. 

Miss Martha E. Graham of the Douglas Aerodynamics Depart- 
ment carried out many special computations and supplied valuable 
criticism. At the Jet Propulsion Laboratory, George Morikawa 
worked on the many script and H. J. Stewart read it critically. 
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INTRODUCTION AND SUMMARY 

The theory of conical flow was initiated by Busemann in Reference 
and developed further in Reference 4 where the foundations of a theory of 
conical wings are laid. However, these references contain only a brief 
outline of the methods and the solutions for only two special cases of 
conical wings. Some powerful methods were introduced but often presente 
with insufficient explanations. The present report* tries to clarify these 
methods and develop them further. It also treats a large number of cases 
not considered by Busemann. The emphasis is on a detailed treatment of 
the general methods and of the basic conical solutions. At the same time 
some new ways of using conical solutions in applied wing theory are 
pointed out. 

Section I is devoted to the foundations of conical theory. Boundary 
conditions are discussed in detail, as well as general methods like the 
oblique transformation. Formulas are derived which will simplify the work 
in the later sections. 

Sections II to VI, inclusive, are devoted to detailed studies of vari- 
ous basic conical solutions. In general, expressions for the entire flow 
field are derived. In the simplest application of conical theory one di- 
vides a given wing of finite chord and span into various regions, in each 
of which some conical solution is valid. The effect of angle of attack is 
also separated from the effect of thickness with the aid of the superposition 

*This report was written under the sponsorship of the Ordnance Depart- 
ment, U.S. Army, Contract No. W04 -2 OO-ORD-4 55 . 
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principle. Basic solutions useful for such purposes are contained in 
Sections II and III. With the aid of these solutions the flow field near 
wings of simple planforms and profile may be determined. This is the way 
conical theory is customarily applied (Cf. Refs. 1, 6, and 7). However, 
the use of conical flow theory goes far beyond this obvious application. 

One further method of application was indicated by Busemann, namely, a 
method of determining flow conditions (downwash and sidewash) behind a 
three-dimensional wing. The method uses wings of constant-lift distribu- 
tion. The solutions for such wings are worked out in Section IV. Refer- 
ence 2 contains some applications to wings of practical interest where the 
method is illustrated in detail. 

Another type of wing, which to the author’s knowledge has not been 
considered before, is the mixed type consisting of a region of constant 
lift adjacent to a flat plate. The solutions are easily obtained with the 
methods developed in the present report. Essentially these are solutions 
of a certain type of interference problem. With their aid the lift problem 
may be solved for a much larger class of planforms than was possible with 
the solutions from Section III. In particular, wings with low aspect ra- 
tio may be treated, as pointed out in Section III-K. 

In Section VI a different type of interference problem is studied. 

The wings considered consist of adjacent flat regions at different angles 
of attack, and have many applications. The application to control surfaces 
has been worked out in Reference 3. 

At the time the essential part of the work in the present report was 
done, conical methods frequently provided the fastest and sometimes the 
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only way of solving many problems in applied wing theory (Cf. Refs. 1, 2, 
and 3). Although other methods have been developed since, it is believed 
that the solutions presented in the present report still have considerable 
practical application. The method of superimposing conical flow fields 
also has the advantage of being very intuitive. Furthermore, the theory 
of conical flow is of considerable theoretical interest. In the Busemann 
theory, the problem of finding the flow around a conical wing is reduced 
to that of finding certain analytic functions of a complex variable. With 
the aid of the powerful methods of the theory of analytic functions, the 
essential structure of the solution may often be seen directly from the 
boundary conditions. In this way a deeper understanding of the nature of 
conical flow fields is provided. It is to be hoped that eventually similar 
results may be obtained for more general flow fields. 

I. GENERAL THEORY OF LINEARIZED SUPERSONIC FLOW * 

In this section some fundamental properties of linearized supersonic 
flow will be studied. In particular, the concepts of conical flow and 
conical wings will be introduced. ‘The main purpose, however, is to furnish 
the necessary basis for the study of specific conical wings in the re- 
mainder of the report. 

A. Fundamental Equations and Boundary Conditions 

In general, it is assumed in this report that the free stream Mach 
number has the value M =s/z. As will be shown in Section I-M, a problem 
at a different supersonic Mach number may be reduced to an equivalent 
♦The nomenclature used in this report is given in Table I. 
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problem at M by a Pr andtl-Glauert transformation. Although this 

reduction is not universally valid, it is admissible in the linearized 
theory for the type of bodies (wings) treated in the present report. The 
free stream velocity will be denoted by w^, taken in the direction X3. 

The fundamental equations for linearized supersonic flow at M =>TT 

are 


8xx xi + 8x 2 x 2 “ 8x3x3 = 0 (g = u, v, w or $ ) 


( 1 . 1 ) 


and the equations of irrotationality interconnecting the various pertur- 
bation velocities and the velocity potential are 

grad $ = (u, v, w) (1.2a) 

curl (u, v, w) = 0 (1.2b) 

Here (u, v, w) denotes a vector (the perturbation velocity) with three 


components u, v, and w. 

If the free stream Mach number is M, Equation ( 1.1 ) will have to be 
replaced by 

8x1X1 + 8 x 2 x 2 - —2 g X3X3 = 0, m = (U,) 

The pressure and the acceleration potential (Cf. Ref. 8) are obtained 

from the perturbation velocity w by 

P " Poo = w c w ( 1 . 3 a) 

acceleration potential = w^ w ( 1 . 3 b) 

When there is no danger of misunderstanding, the component w will be 
referred to as pressure. 

The boundary conditions usually considered are: 
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On the surface of a body in a supersonic stream the 

flow vector (u, v, w^ + w) is tangent to the body. (1.4) 

At a great distance upstream of the body (infinity), 

the velocity vector and the thermodynamic quantities 

have the constant values (0, 0, w^), p^, , etc. (1.5) 

In the problems treated in the present report, condition (1.5) still 
applies. But condition (1.4) will be modified for several reasons: First, 
it presupposes that the shape of the body is given. In some problems, 
however, the pressure will be given instead. In these cases condition 
(1.4) will be replaced by a suitable condition on w. Second, when the 
shape of the given body is of a certain type, condition (1.4) may be sim- 
plified and stated as a condition on v only. Third, in treating conical 
flow some important mathematical transformations of the boundary condi- 
tions will be made which, however, will not change their content or 
introduce any new approximations. In some cases it will also be necessary 
to add some further conditions in order to determine the solutions 
uniquely. 

These modifications of the boundary conditions will be discussed in 
Sections I-B and I-K. 

B. Wings (Planar Systems), Symmetry Properties, and Boundary Conditions 

All specific solutions in Sections II to VI, inclusive, will be for 
conical wings as defined in Section I-K. The fundamental properties dis- 
cussed in Section I-B will apply to conical wings as well as to more 
general wings; in particular, a large class of finite three-dimensional 
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wings as actually used on airplanes or missiles. 

In linearized theory it is often permissible to apply the boundary 
condition ( 1 . 4 ), or some substitute therefor, to a region which is con- 
tained in the XjX 3 -plane. Such a body will be called planar, and in this 
case the X^-plane will be referred to as the plane of the wing. Of 
course the wing is not contained in this plane except for the trivial case 
of a flat plate at zero angle of attack. Prescribing a quantity on the 
wing will then mean prescribing it on the vertical projection of the wing 
on the XjX 3 -plane. 

The simplest case is that of a wing whose shape is prescribed in such 
a way as to satisfy the following condition: 

The tangent plane is everywhere almost horizontal, and 

the wing has some points in common with the horizontal 

XiX 3 -plane. {!’&) 

For such a wing one defines: 

Local angle of attack at the point P = (x lf x z , x 3 ) is 
the angle between the tangent plane at P and the X^X 3 ” 
plane measured in a plane x^ = constant. (1*7) 

For a wing satisfying condition (1.6), condition (1.4) is replaced by the 
following 

If the local angle of attack at P = (xp x 2 , * 3 ) is Ct> 

then v = -w^CC at P = (xi, 0, x 3 ). (I- 4 ') 


page 6 



MCA TN 1685 

JEk 


Progress Report No. 4-3 


Note that here a condition at a point P on the wing is replaced by a 
* condition at its projection P in the plane of the wing. Condition (1.4* ) 

may be expressed as v = known function of xj and x^ in vertical projection 
of wing on -plane. 

When the pressure or the lift is prescribed, condition (1.4 ) will be 
replaced by a condition on w; and in all cases of interest in the present 
report the values of w will be given in the plane of the wing. Then con- 
dition (1.4 ) will be replaced by: 

w = known function of (xj, x^) in a region which is 

the projection of the wing on the X 1 X 3 -plane. ( 1,4 ") 

4 Finally, certain mixed types of wing are of importance. Here v is 

given over one part of the wing and w over a different part. 

Actually, the region in the X 1 X 3 -plane may have a top side and a 
bottom side. Points on the -top side will be considered as projections of 
points on the top surface of the wing. Their X£ -coordinate will be denoted 
by 0+. Similarly X£ = 0- for points on the bottom side. The value of a 
function f at (xj, 0 +, X 3 ) will be considered as the limiting value 

lim f(x 1 , S , x 3 ) 

S 0 

In the same way f(xj, 0-, x^) is defined as 

lim £(xj, -S , X 3 ) 

6->0, S > 0 
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Very often one of the velocity components will have different values 
for = 0+ and = 0-. In particular, such a discontinuity may be pre- 
scribed in the boundary conditions (1.4 1 ) or (1.4 U ). Discontinuities may 
also occur in the -plane off the wing. 

In general, condition (1.5) plus a prescription of the values of v or 
w on the wing will not determine the flow field uniquely. One will need 
conditions on the part of the X 1 X 3 -plane off the wing and sometimes some 
conditions on the behavior at the edges. 

For a study of the first questions some observations are needed about 
symmetry properties of the flow and discontinuities in the X\X$ -plane. * 

Proofs of some of the assertions made later in Section I-B would have to 
be based on the invariance of the Equations (1.1) and (1.2) under reflec- ^ 

tion (Cf. Sections I-C and I-D). 

When the word symmetry is used without specification, itiefers to 
symmetry with respect to the X$Z $ -plane since for the present purpose this 
is the important case. 

A function f(xj, x^, X 3 ) is said to be symmetric, with respect to the 
XjX^ -plane, if it is even in x^; i.e., if 

f(xj, - * 3 ) = f(Xp x 2 , X3) ( 1-8a )’ 

It is antisymmetric, with respect to the X 1 X 3 -plane, if it is odd in x^; 
i.e., if 

£ ( x l» " x 2 * x 3) = -f ( x l» x 2 « * 3 ) (!- 8 b) 

If P = (xj, x 2> X3), then (xj, - x^, X 3 ) is called its reflected point P r . 
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The following simple property is fundamental: 


If f is antisymmetric, it is either zero or dis- 
continuous for X 2 = 0. 


(1.9) 


In proof of this statement, let 


f(x lt 0+, X 3 ) = lim £(x 1( S , X,) = +a 

o-> 0 t 6 ?o 


Then, by condition (1.8b) 

f(x,, 0-, x 3 ) = lim f (x,, - <5, Xo) = -a 
<5-*0 <5> 0 1 

Obviously, a / 0 means a discontinuity of magnitude 2a. With respect to 
symmetry properties there will be two main types of body shapes. 

Symmetrical nonlifting wings . Such wings are bisected by the XjX 3 - 
plane. Whenever a point P is on the upper surface, its reflected point 
P r is on the lower surface. According to conditions (1.4*) and (1.8b), 
a discontinuous antisymmetric distribution of v is prescribed on the wing 
(i.e., on its projection in the X^X^-plane). 

Lifting wings of zero thickness. Here the prescribed distribution of 
v is continuous and symmetric. 

Corresponding to these two types of wings there are two types of flow. 

These two types will be denoted, rather arbitrarily, as symmetric and 
lifting. The properties are summarized in the following tabulation in 
which S - symmetric and A = antisymmetric: 
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Case 

$ u 

V 

w 

body ] 


Symmetric 

S S 

A 

S 

symmetric / 

( 1 . 10 ) 

Lifting 

A A 

S 

A 

zero thickness ) 



The importance of these two types of flow is that any flow may be 
decomposed into a sum of the two types of flow, for may be represented 
as a sum of an odd and an even function, and the behavior of $ determines 
the other quantities uniquely by condition ( 1 . 2 a). 

For this reason all wings considered later in the present report will 
be either symmetrical and nonlifting or lifting and of zero thickness. 

The fundamental principle for determining the boundary conditions of 
the wing is: 

The X 1 X 3 -plane may carry a discontinuity in v or w only 
where there is a wing. ( 1 . 11 ) 

Combining condition (1.11) with conditions (1.9) and (1.10), one 
obtains: 

For a symmetrical wing v = 0 in the X^X-^ -plane off the 

wing. ( 1 . 12 a) 

For a lifting wing of zero thickness w = 0 in the 

X 1 X 3 -plane off the wing. ( 1 . 12 b) 

Note that a discontinuity in u is possible off the wing. This dis- 
continuity is nothing but a vortex sheet. 

Reference 2 contains a detailed discussion of some of the principles 
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which have just been mentioned. They are actually of great importance for 
supersonic wing theory for problems of pressure distribution in an induced 
flow field behind a finite wing. They will also be used consistently in 
the present report both in the derivation of the basic solutions and in 
Sections IV-B, V-E, and VI-E. 

Conditions (1.12) and (1.4) prescribe either v or w over the entire 
XjX^ -plane. There may then occur two types of boundary value problems 
which are essentially different from a mathematical point of view. Take 
a symmetrical wing with given shape, say, of constant angle of attack. 

Then v = + v Q on the wing and v = 0 off the wing in the X^X^ -plane. Hence 
v alone is prescribed. This problem is essentially easy. Although 
Busemann’s conical methods may furnish rapid and elegant solutions (Cf. 
Section II), there are other straightforward solutions, e.g., by a source 
sheet which may be determined directly from the shape of the body. Simi- 
larly in the lifting case if the lift is prescribed, w = + f(xpX 3 ) on the 
wing where f is a known function and w = 0 off the wing. Again this is a 
straightforward problem. The solution may, e.g., be had from the accelera- 
tion potential (Cf. Ref. 8) which is directly given from the lift distri- 
bution. However, conical methods may again be used to great advantage 
when the planform is conical (Cf. Section IV). Finally it should be noted 
that since v and w satisfy the same equation, any solution for a symmetri- 
cal wing of prescribed shape yields a solution for a wing of zero thickness 
and prescribed lift distribution if v is replaced by w. 

However, if a lifting wing of prescribed shape is given, the problem 
is much more difficult. On the wing v is prescribed; off the wing w is 
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known (= 0). Most theories of a lifting surface have led to an integral 
equation which is very easy to set up but for whose solution there exist 
no straightforward methods (Cf. Refs. 9 and 10). In the solution of such 
a problem Busemann's methods show their greatest power. The lift problem 
is reduced to one in classical potential theory which may be solved by the 
powerful tools of theory of complex variables. Ways by which the mixed 
boundary value problem may be reduced to well-known problems will be dis- 
cussed in Section I-K. Of course, Busemann's method applies directly to 
conical wings only, but one may obtain a very large class of finite wings 
by superposition of conical wings (Cf. Refs. 1, 2, and 3 or Section V-E). v 

Finally for a mixed type of wing (Cf. Section V) the situation may be 
even more complicated. For the lifting case v is prescribed on part of * 

the wing, w is prescribed on some other part, and w = 0 off the wing. 

Again Busemann's concept will provide easy solutions. 

In some applications the addition of condition (1.12) is not sufficient 
to make the solution unique, but it is necessary to prescribe the behavior 
of the solution at the edges. First, some definitions are needed. An 
edge is called supersonic if its sweepback angle is less than the comple- 
ment of the Mach angle; i.e., if the projection of the free stream velocity 
Wqo on a plane normal to the edge is still supersonic. If this projection 
is subsonic or zero, the edge is said to be subsonic. Edges may also be 
classified as leading and trailing. If a line in the X 3 -direction leaves 
the wing when crossing an edge, the edge is said to be trailing; in the op- 
posite case it is leading. In the intermediate case when the edge, or 
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rather its projection on the X^X^ -plane, is parallel to the x 3 -direction, 
it is called an unyawed side edge. The essential assumption about the 


edges is: 


The behavior of the flow at a subsonic edge is 
qualitatively the same as in incompressible flow. (1.13) 


This rather vague principle will new be amplified. Consider a lead- 
ing edge. If it is supersonic, conditions at the edge are easily deter- 
mined from the geometry of the wing and of the Mach waves (Cf. Section 
I-K). However, at a subsonic edge the flow singularity is of the type 
encountered in subsonic flow. In the lifting case the perturbation ve- 
locity becomes infinite as the distance to edge to the power (-1/2). In the 
symmetrical case the edge represents a logarithmic singularity (e.g., Cf. 
Sections II and III; and Ref. 11 for a general discussion of the principles). 

Now consider a trailing edge in the lifting case. As the flow crosses 
the edge, w has to become zero. This is the principle of cancellation of 
lift behind the wing (Cf. condition 1.12b). If the edge is supersonic, 
the lift may be cancelled discontinously by appropriate Mach waves at the 
edge (Cf. Ref. 2 for a discussion of the implications of this possibility). 

At a subsonic trailing edge, however, a finite discontinuity is not possible. 
There are two mathematical possibilities: Either the magnitude of w de- 
creases continuously on the wing to zero, or there is a discontinuous shift 
from 00 to 0. The latter case would simply be a reversal of what happens 
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at the leading edge. In accordance with condition (1.13) it will be as- 
sumed that the first possibility corresponds to the actual behavior of the 
fluid (Kutta condition). The following two implications of condition 
(1.13) will actually be of decisive use in solving special problems in 
this report: 

On a flat lifting wing, | w| decreases continuously to zero 

at a subsonic trailing edge and increases to infinity 

as the distance to the (-1/2) power at a subsonic leading 

edge (Cf. Section III). (1 . 13 f ) 

Near a subsonic edge (leading or trailing) of a flat 

symmetrical wing, |w| has a logarithmic singularity 

(Cf. Sections II-A and II-C) . (1.13”) 

C . Invariance of Fundamental Equations Under Transformations 

The concept of the invariance (or, more generally, of the transforma- 
tion) of an equation under a mapping or transformation of variables is 
very important in physics and mathematics in general. In Sections I-C to 
I-G inclusive, the behavior of the linearized supersonic equations under 
certain mappings will be investigated. This study will furnish an under- 
standing of the symmetry properties discussed in Section I-B , a method for 
treating yawed wings, and finally a basis for the theory of conical flow. 
Similar concepts will be used later in treating the Prandtl-Glauert trans- 
formation (Cf. Section I-M) . On the whole there are two different types 
of results from an investigation of transformation properties: (a) certain 
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symmetry properties of solutions and (b) methods for generating new solu- 
tions from a given one. 

The invariance of an equation under a transformation may be explained 
in two different ways: 

In the first interpretation and xf are considered to be different 
systems of coordinates for the same space. If the equation does not 
change its form when the coordinates x* are introduced instead of x. ( then 
it is said to be invariant under the change of variables x^— ^x* . Equa- 
tion (1.1) would have to be 

Kx l* x l* x 2 * x 2* x 3* x 3* 

in order to remain invariant. 

In the second interpretation x* is considered to be a mapping 

from a space S to a space S* whereby the point P = (xj, X£, X3) is mapped 
on the point P* = (x^*, x^*, xy*). This correspondence between points in 
S and S* induces a correspondence between functions defined in the two 
spaces. To a function f in S there corresponds a function f* in S* defined 
by the rule that f* takes the same value at P* as f at P. Similarly, a 
relation or equation in S will correspond to a transformed relation or 
equation in S*. Examples of such correspondence are given in Sections I-F 
and I-M. If the transformed equation is invariant, i.e., the same as the 
original equation, one obtains a new solution f* of this equation for each 
given solution f. Thus the mapping is used for generating new solutions 
(this is the way the oblique transformation is going to be used). Even the 
case when the equation is not invariant is important. Here the mapping is 
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used to generate solutions of the transformed equation (Cf. Section I-M 
for discussion of Prandtl-Glauert transformation). 

Of these two interpretations of invariance, the second, though less 
familiar, is definitely preferable in order to put the discussion on a 
firm basis. The abstract presentation just presented will be illustrated 
by examples in Section I-D to G. 

The following four transformations, which are of interest for the 
present report, will not be discussed: 

1. Reflection in the plane of the wing. 

2. Rotation around the X 3 ~axis. 

3. Lorentz transformation in the plane of the wing (oblique 
transformation) . 

4. Uniform expansion around the origin. 

D. Reflection in the Plane of the Wing 

By reflection in the plane of the wing is meant the transformation 

— >x*, where x^* = +x^, ~ ~ x 2 » x 3 * ~ x 3 (1-14) 

Equation (1.1) is obviously invariant under this transformation. Let P 3 ^ 
(xj.-x^, x-j) be the reflected point of P = (x^ x^, X3). Let $ (xj, x^, x^) 
be a solution of Equation (1.1) and define §>* as the function whose value 
at P* is that of $ at P. In other words, (x 1 rx 2 ,x 3 ) = $ (xj.x 2 .x 3 ). 

In accordance with the general reasoning in Section I-C, is also a 
solution of Equation (1.1) because of the invariance under reflection. 
Applying Equation (1.2), one sees that if the perturbation velocities at P 
are u, v, w, then at P* the corresponding quantities u* = u,-v* = v, w* = 
Page 16 
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Thus Equation (1.2) is invariant under transformation (1.14) only if v is 
replaced by v* = -v in the mapping. In this way a solution $ has generated 
another solution $*. If ^is even in x^, then the two solutions are iden- 
tical; i.e. $ (P*) = & *(P*), u(P*) = u*(P*), etc., and there results the 
case of a symmetrical flow discussed in Section I-B. The important thing 
is that, conversely, if the boundary conditions are consistent with sym- 
metric flow (e.g., if the shape of a symmetric body is given), it follows 
from the invariance discussed that the solution also is a symmetric flow. 
Actually, this conclusion is justified only if the boundary conditions 
determine the solution uniquely. The logic of this proof of symmetry will 
be discussed in greater detail in connection with conical symmetry (Section 
I-G). 

Similarly, since Equation (1.1) is homogeneous, - <J)* is a solution if 
<£)* is. Hence the solution ^(xj, X£, x^) generates the solution - $ (x^, -x 2 ,X 3 ). 
If the original is odd in x^, the two solutions are identical. This is 
the lifting case of Section I-B. 

Of course, one may also study reflection in the -plane and the 

corresponding symmetry with respect to this reflection. This case is, 
however, much more obvious. The main result for symmetries with respect 
to this plane is that, if the body is symmetrical, then v and w are sym- 
metrical and u antisymmetric. In particular 9w/3xj = 0 in this plane. 

E. Rotation Around the X 3 -Axis 

Obviously both Equations (1.1) and (1.2) are invariant under rotation 

around the X 3 -axis. This is the basis for studying bodies with rotational 
symmetry. It is mentioned here only as a preparation for introducing the 
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oblique transformation. The significant fact is that Equation (1.1) is 
not invariant under a rotation about any other axis . 1 Thus the x^-axis is 
a preferred direction, whereas the xj- and x^-axes are equivalent. In 
incompressible flow, on the other hand, the equation corresponding to 
Equation (1.1) is obtained from Equation (1.1) by changing the minus sign 
ahead of the last term into a plus sign. The resulting Laplace equation 
is invariant under all rotations. One might then ask whether there are 
transformations which are analogous to general rotations but which leave 
Equation (1.1) invariant. 

Now an ordinary rotation around a point Q is a transformation P-»P* 
which leaves the distance to Q unchanged: | QPl =|QP*| . The answer to 
the question above will be found by introducing a new metric based on the 
concept of hyperbolic distance. A hyperbolic rotation around Q, also called 
a Lorentz transformation, will be a transformation which leaves the hyper- 
bolic distance to Q unchanged. These concepts will next be defined and 
discussed. 

F. Hyperbolic Distance, Lorentz Transformation, Oblique Transformation, 
and Y aw 

The following definition of hyperbolic distance between the points 
P = (xp x^. x^) and P= (xp x^, X 3 ) will prove convenient for dealing 
with supersonic flow when the Mach number is equal to {~Z and the free 
stream is parallel to the X 3 -axis. x 

_ r _ 2 _ 2 _ 2 f 

Hyperbolic distance PP = I (xj-x j) - (x^ - x^) - (xj - Xj) J (1.15) 

Note that the Mach cone from a point P may be defined as the locus of all 
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points whose hyperbolic distance to P is equal to zero. 

A hyperbolic rotation around Q, or a Lorentz transformation, is de- 
fined as a transformation P-»P* such that the hyperbolic distance QP = 
hyperbolic distance QP*. It follows that Q is mapped on itself and that 
the Mach cone from Q is invariant in the sense that any point on this Mach 
cone is transformed into some other point on the same cone. 

The following classical theorem will be the basis for the discussion 
in the remainder of Section I-F. It can easily be proved directly for the 
special cases which will actually be used. 

Equation (1.1) is invariant under any Lorentz 

transformation. (1,16) 

To obtain a more explicit description of these transformations, con- 
sider first the case where X 3 * = x^. As is easily seen, invariance of 
hyperbolic distance is then equivalent to invariance of ordinary Euclidean 
distance. Thus such a transformation is simply an ordinary rotation 
around the x^-axis, and theorem (1.16) reduces to the invariance discussed 
in Sectiop I-E. 

However, when X 3 ^ X 3 *, the hyperbolic rotation differs from the or- 
dinary rotation. One will need only the case of a hyperbolic rotation 
around the origin which leaves the plane of the wing invariant. Such a 
special Lorentz transformation will be called an oblique transformation 
and its explicit expression is: 
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The fact that this transformation leaves both the hyperbolic distance 
and Equation (1.1) invariant is easily checked. Theorem (1.16) depends, of 
course, on the fact that Equation (1.1) is a special case of the wave 
equation. The general equation for a disturbance (e.g. optic or acoustic,) 
spreading in (x,y, z)-space with velocity c is: 

8 xx * Syy + Szz ” c 2 ®tt ~ ^ (^-*^) 

This equation reduces to Equation (1.1) in the obvious way. In other 
words, if X 3 is thought of as denoting time, Equation (1.1) is the equa- 
tion for a disturbance spreading in the X^X^-plane with unit velocity. 

This analogy was first noted and exploited by von Karman (Cf. Ref. 12). 

That Equation (1.18) is invariant under a four -dimensional Lorentz 
transformation is a well-known fact from mathematical physics and is funda- 
mental in wave theory and theory of relativity. In aerodynamics it has 
been used by Kuesner and others (Cf. Ref. 13) for the study of nonstation- 
ary processes like fields due to moving sources. When Equation (1.18) is 
specialized to Equation (1.1), the general theory of Lorentz transformation 
yields the assertion (1.16), in particular the invariance of Equation (1.1) 
under Equation (1.17). Although this is an obvious consequence of von 
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Karman’s analogy, to the author’s knowledge it has never been utilized in 
aerodynamics until recently when Jones (Cf. Ref. 14) used the oblique 
transformation (1.17) to treat a special case of a yawed wing. The idea 
of treating yawed wings by this method has proved quite fruitful; it will 
be discussed further and applied ofteh in the present report. 

Some formal properties of Equation (1.17) will now be developed. 
Since = x^*, the oblique transformation affects only the plane of the 
wing and may be represented by a two-dimensional matrix 



The Jacobian of this transformation is the determinant 
of M and is unity. Hence the oblique transformation 
is area-preserving. (1-19) 


That this statement is of importance for evaluating lift of a yawed wing 
will be clear from the discussion later in Section I-F. 

The inverse transformation of Equation (1.17) is given by 
xi* - axs* 


x 2 = x 2* 


( 1 . 20 ) 


-ax]* + X3* 
n/ 1 - a 2 
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These equations mean that the inverse of M in Equation (1.17 ) is given by 


1 / 1 


(1.20') 


V 1 - V -a 1 / 

Now let $(x r x 2 . x 3 ) be a solution of Equation (1.1). Then in accord 
ance with the general principle discussed in Section I-C, one may generate 
another solution from $. Let P* be obtained from P by the transformation 
(1.17). Then put 


$*(P) = $(P*) 


(1.21a) 




Using Equations (1.17), Equation (1.21a) may be written 


_ / xi + axi axi + x-A 

$. (W3 ) . 7 =) <w«» 


If u* is the sidewash belonging to the generated solution 




3x i*\ 


9x i / ' vT^ 


= [u(P*)+a w (P*j) 


Carrying out the same computation for the other velocity components 


yields 


u*(P) = . 1 - ju(P*) + a w(P*)J 

V 1 - a‘ L 
v*(P) = v(P*) 

v*(P) = 1 [a u(P*) + w(P*)] J 

n/ 1 - a z 



(1.21c) 


^Generally in this report u*, v*, and w* will denote the harmonic conjugates 
of u, v, and w (Cf. Section I-J). The notation in Equations (1.21) and (1.22) will 
not be used later on. 
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Thus equation (1.2) is not invariant under Equation (1.17). As a 
matter of fact, Equation (1.2) is instead invariant under ordinary rota- 
tions. In general, w*(P) = w(F) contradicts Equation (1.21c). However, 
since w also obeys Equation (1.1), one may define a new solution w* by 

w*(P) = w(P*) (1.22) 

But if u* and v* are then computed from Equations (1.2) and (1.22), they 
will in general not obey Equation (1.21). 

In the applications both schemes (1.21) and (1.22) will be used de- 
pending on the problem, the choice of scheme having to do with the bounda- 
ry conditions at the edges. The discussion of this and some other 
important features of the oblique transformation will be presented in 
Section I-L. 

G. Uniform Expansion. Conical Flow 

By a uniform expansion around the origin is meant the transformation 
x i — *x.* = CX 4 , c > 0 (1.23) 

Evidently Equations (1.1) and (1.2b) are invariant under this trans- 
formation, but Equation (1.2a) only if £ is replaced by c <J>. 

Invariance under reflection was related to certain symmetry properties 
as discussed in Section I-C. Similarly the invariance of the equations 
under Equation (1.23) will imply a certain conical symmetry of the flow 
wherever the boundary conditions themselves are invariant. 

A flow field is said to be conical with respect to a point P (apex) 
if u, v, w, and thermodynamic quantities like pressure, density, etc. are 
constant along any half-line (ray) issuing from P. Unless otherwise 
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specified, the apex will always be assumed to be the origin of the system 
of coordinates. Then the definition of conical flow above may be re- 
phrased: A flow field is conical if it is invariant under Equation (1.23). 
Finally one may say that a flow field is conical if u, v, etc. depend only 
on the ratios X 1 /X 3 , x^ A 3 - 

A conical body with apex at the origin is formed when a half-infinite 
straight line with a fixed end point at the origin moves along some space 
curve. The simplest example is of course the circular cone, but an infi- 
nite flat triangle with apex at the origin is also a conical body. 

Boundary conditions are said to be conical whenever they are invari- 
ant under Equation (1.23). Thus the same values of perturbation velocities 
or whatever is given must be prescribed for the point (xj, X£, X3), and any 
point (cx^, cx^, CX 3 ) where c > 0. The simplest example of such boundary 
conditions occurs when the shape of a conical body is given. The other 
important example for the present report occurs when the pressure is pre- 
scribed in the plane of the wing in such a way that it is constant along 
each ray from the origin. 

The basis for the theory of conical flow is the following assertion 
about conical symmetry; 

A flow field is conical if the boundary conditions are 

conical. (1.24) 

This statement is actually valid for the full nonlinear Euler equations 
but fails when viscosity is introduced. Because of its fundamental impor- 
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tance the proof will be considered in some detail even though it may be 
claimed to be rather obvious. It may be proved in two different ways 
which are at least seemingly different: first, a physical argument based 
on dimensional analysis and, second, a more mathematical analysis depend- 
ent on the principles of invariance discussed in Section I-C. 

The first argument runs as follows: In a nonviscous flow no combina- 
tion of the quantities given at infinity (w^, P<^ etc.) has the dimension 
length. Furthermore, a conical body or any conical boundary condition is 
described without specifying any length. Thus the problem contains no 
characteristic length, and the only way to form nondimensional combinations 
involving the coordinates is to divide these by each other; for example, 
u/w 00 has to be a function of x^ /x^ , X 2 /x^ , and some nondimensional combi- 
nation of the values at infinity. 

The second method of proving Equation (1.24) may seem more involved 
but will be given as a concrete example of the principles discussed in 
Section I-C. It should be emphasized that this method works in many cases 
where dimensional analysis fails (e.g., for the Prandtl-Glauert transfor- 
mation). Consider a stationary supersonic flow around a body B (not 
necessarily conical) in a space S. Assume that the solution for this flow 
has been obtained. Now construct a new solution in another space S* by 
the following rule: To each point P in S let a point P* in S* correspond 

where P* is obtained from P by the transformation (1.23). The body B in 
S will then correspond to an expanded body B* in S*. Define a flow field 
in S* by letting the velocity at P* be the the same vector as that at P and by 
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letting conditions at infinity be the same in both spaces. 

The flow thus defined is the correct solution for the flow around B* 
in S* with the given conditions at infinity; therefore Equations (1.1) and 
(1.2b) and the boundary condition (1.4a) are fulfilled. The first part of 
this assertion follows from the invariance of the equations under the map- 
ping rule (1.23). The second part follows from the fact that the normal 
at F* on B* has the same direction as the normal at P on B. 

Now let the body be conical. Then B and B* are the same (or rather 
congruent) bodies. Hence assuming that the problem (equations and boundary 
conditions) has a unique solution, the flow field in S must be exactly the 
same as in S*. But from the way the flow field in S* was constructed, it 
follows that flow conditions at the point ( cxj, cx^, cx^) in S must be the 
same as those at (x : , x 2 , x 3 ). Since this statement is true for any posi- 
tive c, the flow in S must be conical. 

Thus assertion (1.24) is proved for the case when the shape of a 
conical body is prescribed. The proof for any type of conical boundary 
conditions is of course similar. 

As is well known (Cf. Ref. 4), the above considerations apply also to 
$ /x_j (but not to $ itself). However, very little use will be made of this 
fact in the present report. 

H. Reduction s of Equations for Linearized Supersonic Conical Flow. 

Tschaplygin Transformation 

Although the reasoning in Section I-G was valid for the nonlinear 
Euler equations, the discussion will now be restricted to the linearized 
equations. As pointed out in Section I-G, in conical flow the number of 
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independent space variables may be reduced from 3 to 2 . This reduction 
will now be carried through explicitly for Equation ( 1 . 1 ). 

Consider a solution g(xj, x^, x^) of Equation (1.1) where g is either 
u, v, or w (but not $). To say that g is invariant under Equation (1.23) 
is the same as saying that g is homogeneous of order zero in xj, x^, x^. 
Introduce cylindrical coordinates (r, 0 , X 3 ) by 


0 = arc tan x 2 
X 1 


(1.25a) 


r = V X 1 Z + X 2 Z (l-«b) 

Coordinates 0 and r are homogeneous of order 0 and 1, respectively, in x^ 
and x^. By repeated application of Euler’s theorem for homogeneous functions 

7 7 

*3 e *3 = ' lTl X i 8 *i = ‘ ti (x i g 9' e xi + *i» r «- XI ) 

= - 0 -g e - r-g r = - r-g r 

Since x-^g is also homogeneous of order 0, by the same principle 


X 3[ x 38x 3 1x 3 = - r (' r g r ). 


Thus 


x 3 g *3x 3 => I 3[ x 3 8 x s ] x ‘ X 3«* 3 = r < r 8r> r + r ' 


Sr 


(1.26) 


It is also known that the two-dimensional Eaplacian may be expressed 
in polar coordinates as 


1 , 1 

g xixi + §x 2 x 2 = r -(r-g r )+72 g ee 


(1.27) 
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Finally it is observed that a conical flow is completely known if it is 
known on the plane x^ = 1. Hence a reduction of variables may be effected 
by introducing Equations (1.26) and (1.27) into Equation (1.1) and then 
putting Xj = 1. The result is 

( 7- ‘ r ) ( r ‘ g r) r ’ r ’ 8r + 7 8 ee = 0 U- 28 ) 

This is thus the equation for linearized conical flow at the station 
x^ = 1. In this plane the Mach cone from the origin is represented by the 
unit circle (Mach circle). Parallel to the Xj-and X^-axes in the plane 
x^ = 0 , the Xj- and X 2 ~axes may be introduced into the plane X 3 = 1 with 
the same coordinates x^ and x^. In this plane Q and r as introduced by 
Equation (1.25) become polar coordinates. The X^-axis goes through the 
origin of the plane. 

By applying the usual formal criteria (Cf. Ref. 15), one sees that 
Equation (1.28) is elliptic inside the unit circle and hyperbolic outside. 

An intuitive reason for this fact can be found by studying the domains of 
dependence of the original Equation (1.1). An interesting discussion of 
this question is given in Reference 4 . 

Now the question arises: Since Equation (1.28) is elliptic inside the 
unit circle, does there exist a transformation of coordinates which reduces 
it to the simplest of all elliptic equations, namely, the Laplace equation? 
This question was first answered in the affirmative by Tschaplygin, who 
studied a similar equation in a different connection. A function A = Ts(a) 
for any real number between -1 and 1 may be defined by 
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Ts(a) = 


JT- a 2 


(1.-2 9) 


(1.30) 


This function will be called the Tschaplygin transformation of a. The 
change of coordinates introduced by Tschaplygin is then 

0 unchanged ) 
r-*R = Ts(r) J 

It should be emphasized that there is a straightforward way of finding 
this transformation; namely, it is natural, first of all, to leave the 
angle 0 unchanged since both Equation (1.2 8 ) and the Laplace equation are 

invariant under a rotation around the origin. Let R = Ts(r) be unknown. 

[ R* 1 7 

8 RR* rI + R iS R J> where P r irne denotes dif- 
ferentiation with respect to r. Equation (1.28) becomes 


R *(i " r 2 )( 8 RR R + jjT g R ) +(7 - 2r)R‘ g R + 7 


rZ~ 8qq - 0 


( 1 . 28 *) 


In order to be the Laplace equation, formula (1.2 8 ) has to be identically 
equal to 

1 1 

SRR + R 8 R + r2 g 00 = 0 (1.31) 

Equations (1.28 ) and (1.31) are the same if the following two equations 
are true for any r: 

(R')V(l - r 2 ) = R 2 (1.32a) 

r 2 (l - r 2 )R" + r 2 ( r “ 2r)R' = R (1.32b) 

The problem is solvable only if these equations are consistent. This is 
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the case since Equation (1.32b) follows from Equation (1.32a) directly by 
differentiation. And Equation (1.32a) can be written 

R >' ■ 

which immediately integrates to Equation (1.3 0) if the solution is normal- 
ized by requiring Ts(l) = 1. 

Following Busemann, the Cartesian coordinates (x, y) are introduced 
corresponding to (R,0) and also complex coordinates by 


r ie 

S> = xj + ix^ = r • e 

(1.33a) 

... _ i0 

^ - x + iy = R ‘ e 

(1.33b) 


The plane * 3 = 1 with its natural coordinates will be referred to as 
the -plane; and the corresponding plane after the Tschaplygin transforma- 
tion, as the e -plane. The Tschaplygin transformation maps each point 
r . e i0 in the -plane on a point e = R • e in the £ -plane where 
R = Ts(r). Also £ will be denoted by Ts( £ ). This transformation leaves 
the unit circle and the origin invariant point by point and also transforms 
the rays from the origin into themselves, although not point by point; 
e.g., if ^ is real and positive, Ts( £, ) is real and positive but different 
from £, except when £, =1. 

Busemann (Cf. Ref. 4) gives a geometrical interpretation of the 
Tschaplygin transformation and also points out that it may be written 

(L34) 
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In the applications one is often interested especially in flow condi- 
tions in the plane of the wing. Then it is convenient to introduce (Cf. 

Fig. 1) 


t = 


X 1 


and 


T = arc tan t 


~ * (1.35) 

Either one of the coordinates t or f specifies a conical ray in the plane 
of the wing uniquely; e.g., pressure distributions on such wings are func- 
tions of t (or 't) only. In the ^ -plane one has the simple relations 


is on 


the Xj-axis: 


t = r • cos 0 = Xp r = 1 1 | 

The Tschaplygin transformation of t is denoted by T 

1 - J 1 - .t 2 


T = 


(1.36a) 


(1.29 ) 


Then 


T = R * cos 0 - x, 


(1.36b) 



FIGURE 1. COORDINATES ON THE PLANE OF THE WING 
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I. Algebraic Properties of the Tschaplygin Transformation 

Let A be Ts(a). In the application a will be r, t, or tangent of the 
angle of some side edge. Then Equation (1.29) implies the following use- 
ful relations between a and A: 


“2 

1 + A* 

(from Eq. 1.34) 

(1.37a) 

1 2 
A + A = a 

(from Eq. 1.37a) 

(1.37b) 

(1 ± A ) 2 = 1 + i 


(1.37c) 

2A a 


A 2 + ‘ - 2 ( 2 - 


(1.3 7d) 

A 2 a 2 


1 1 + Vl - a 2 


(1.3 7e) 

A a 


2 _ 1 - n/ 1 - a 2 


(1.37f) 

A _ i* fi ~r? 


‘‘A 2 ,J 1 - a 2 

1 + A 2 


(1.37g) 

A a 


(1.37h) 

1 - A 2 2 V 1 - a2 


1 + A / 2(1 £ a) 


(1.371) 

f a ^ a 


Also it may be seen from the binomial theorem that 


A = |a for a<< 1 


(1.38) 

Relations Between the Velocity Components 

in the £ -Plane 



In Section I-H the following basic result taken from Busemann was 


proved: 
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In conical flow, the velocity components u, v, w, 

satisfy Laplace equations when expressed as 

functions of the coordinates x and y defined by 

Equation (1.33b). (1.39) 

This statement means that there must exist complex -valued functions 
U, V, W such that U, V, W are analytic functions of e 

U = u + iu* 


V = v + iv* 


(1.40) 


W = w + iw* 


where the real-valued functions u*, v*, w* are the harmonic conjugates 
(determined only within a constant) of the potential functions u, v, and 
w. From now on, an asterisk on the symbol for a perturbation velocity 
will always denote the harmonic conjugate. The function U, V, W will be 
called the complex velocity components. 

Since u, v, w are interconnected by Equation (1.2), there must exist 
some analogous relations for U, V, and W. These relations should follow 
from Equations (1.1) and (1.2) and the fact that the flow is conical. 
Busemann (Cf. Ref. 4) states the following result: 


dW 


d(u + iv) = -|( — + tdW) 


(1.41) 


This equation may be proved as follows: The right-hand side may be 
written as 


|dw-idw* ig 

+ Re (dw+idw*) 


Re 


-ie 


•* 4 


= e“T(R + ^)dw+i(R - 'A' - 


i0/dw .R 1 , / x 

= -e ( — -i^r^z dw *) ( a ) 
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Use was here made of Equation (1.37b) and the formula obtained from Equa- 
tion (1.37b) by taking the derivative with respect to r: 


r 2 R ' R 


(b) 


Thus Equation (1.41) is equivalent to the two relations (c) and (d): 


9(u + iv) 
9R 


i 0 / 1 0w .1 1 3 w 

_e [7 m+'Tz ^ 


R 3 0 


(c) 


or, multiplying by R , 


3 (u + iv) 


3 r 


ie 1 /?w . 1 

r (57 


(=') 


and 


9 (uji iv) 
9 9 


= -e 


ie 


1 22L -i R 1 m 9 w il 

r 3 0 R‘ r 2 (R 3 R ’] 


(d) 


or, using relation (b). 


9 (u + i^ _ 


00 


iO f 1 

■ b 


9 w _ 2\ 2 w 

r JT ^ l - r) T~r 


(d* ) 


In order to prove relation (c 1 ), consider Equations (1.1) and (1.2) 
and homogeneity 


9 w 

9 u 

9 u 


" 9 x 3 

= -r 0 r 

9 w 

av 

0 v 

9x 2 

-0x 3 

“~ r a r 


Multiply the second equation by i and add 

grad w = -r 


which is relation (c*) 


9 (u + iv) 
0 r 
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For proof of relation (d # ), consider Equations (1.1) and (1.2) 


9u 9 v 9 w 9 w 
2 x^ + 3 ~ 9 x 3 " r 9 r 

3u y v 
9x 2 = 9x x 

Making use of these equations to transform the operator ^r- 

oG 


1 9u 
r 


- sin 8 £iL_ + cos 0 
d * J 



I 'll 

r 30 


= cos 0 


3 v 

3x, 


+ sin 


0£v_ - sin 0*2" 


ax- 


9x : 


= 2 ^ - sin 0 ISL 
dr 9 x 3 


- sin 0 ft . X _ + cos 0 

dx± 


3 x 


(cos 0 ^7-+ sin 0 7777-) + cos 0 ^ w 


3 x, 


3 x 


3 x. 


3 u 3 w 

— + cos 0 

3 r 3 x. 


1 _ 3(u + iv) _ g (v - iu) 3 w 
r 3 0“ 3 r + 3x 3 


(i cos 0 - sin 0) 


. 3 (u + iv) . i 0 

= -l — 1 L + le 

9 r 


9 w 



= -l 


9(u + iv) 
d r 



Now use relation (c'): 
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9 (u + iv ) 


= - n 


[ - e ie I (1?L +i 1 23 , e i9 r . _9w] 
L r '9r +1 r 9 0 + e r gr J 

= -e i0 f-i(l^ + iii^) + ir 2 il] 

L ' 2 r r 96 ' 3r J 

= - e ie r 2 ) lJ?l 

U 3 6 V 3r J 


which is relation (d 1 ). 

It is sometimes convenient to transform Equation (1.41) to the form 
used in Reference 16. 

Taking the complex conjugate of each side of Equation (1.41) yields 

, /dW - 

d(u - iv) = - 2 ( e + * dw ) 


Adding this to Equation (1.41) yields the relation 

2du 


uo luc i ciauuu 

= -if (fc+ i- ) dW + (e + i- ) dw] 


But 


2du = dU + dU 


Hence 



Similarly 


2idv = - 



((- ~ ) dw 


- (6 



and 


2dv = dV + dV 


Hence it has been proved that 

dU 

dt 
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M 


•f 


and 


dV 


1. dW 


de =+ ~ t) de 


(1.42b) 


Relations (1.42) will be called the compatibility relations. Three 

analytic functions U, V, and W give the three velocity components of a 

conical flow only if they satisfy these relations which take the place of 

the equations of continuity and ir rotationality. When one of the three 

functions is given, relations (1.42) determine the other two only within a 

constant. The imaginary part of this constant is entirely arbitrary. The 

real part may be determined by the fact that u = v = w= 0in the region 

of undisturbed flow and by the rules for the change of these quantities 

through a Mach wave or Mach cone (Cf. discussion following Eq. 1.44). 

A very important consequence of relation (1.42b) was pointed out by 

Busemann (Cf. Ref. 4): Put (^?L) = c. If W is regular near the origin, 

d e 

^=0 

then dV/d6 is there of the form 


Hence 


dV ic_ 

de 2 e 


function regular at e = 0 


V =.y- in fc + function regular at e = 0 (1.42c) 

Lr 

Similarly U = - — An e. + function regular at <£ = 0 (1.42d) 

Thus even if W is regular at the origin, U and V have logarithmic singu- 
larities there if c^O. In particular, if c is real, i.e., (dw*)^_Q = 0, 
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then the upwash has a finite jump at the origin; v(P) increases suddenly 
by (cTT)/2 if P passes through the origin on the X-axis in the negative 
direction. 

This observation is of great importance in finding functions which 
satisfy the correct boundary ‘Conditions on the wing. It will be further 
discussed in Sections I-K and I-L, and examples of its application will be 
given in Sections II, V, and VI. 

In the plane of the wing e = T. Hence, by Equation (1.37a) 

i U + L)_i±ii = 1 

2 ' ' 2T t 


Thus in the plane of the wing 


du 

dt 


JL 

t dt * 


t = tan T= x 


1 


(1.43a) 


Similarly on the imaginary axis £ = iy and 



1 1 + v Z 

2 y 


1_ 

x 2 


Hence in the X^X^ -plane 


dv JL dw (1.43b) 

dx 2~ x 2 dx 2 

From these formulas it follows easily that if w is constant in a 

cdnical region in the plane of the wing, tl is also constant there. More 

generally, even when the flow is not conical, it follows directly from the 

equation of ir rotationality that U is constant along a (finite or infinite) 

line segment parallel to the X^-axis as long as2w/2X } = 0 on the segment. 

However, if w has a discontinuity on such a line, u might also change dis- 

continuously. If the segment extends upstream into the undisturbed flow, 
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the constant value is of course 0. Behind a trailing edge the constant 
might be different from zero, which for the lifting case means that there 
is a vortex line of constant strength. Examples of regions of constant w 
are, in the lifting case, the part of the plane wing not occupied by the 
wing and regions on the wings studied in Sections V and VI (Cf. Chap. 8C 
of Reference 2). 

K. Conical Wings: Classification. Boundary Conditions in Physical Space 

and in the £ -Plane 

If a conical body is at the same time a wing as described in Section 
I-B, it will be called a conical wing. These are the only conical bodies 
studied in the present report. It will always be assumed that they are 
planar systems, i.e., that the boundary conditions may be applied to the 
plane of the wing. Of course one might also consider cases like wings 
with dihedral, but they will not be treated in this report. Such cases 
may easily be reduced to planar systems (Cf. Ref. 10). 

Conical wings will be classified in various ways; for example, in 
Section I-B they were classified with respect to symmetry conditions of the 
flow or the nature of the quantities prescribed on the boundary. When the 
shape is given, it will generally be assumed that the wing is flat; i.e., 
v = constant. There are two such cases: (1) In Section II flat symmetric- 
cal wings will be treated. These have wedge-shaped profiles of constant 
half angle A , the word profile meaning the intersection of the wing and a 
plane x\ = constant. (2) In Section III the flat lifting wings of zero 
thickness are treated; the profiles are here straight lines of constant 

angle of attack & . After the flat wings, the wings with given pressure 
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distribution will be considered. Since the lifting case is of most inter- 
est, in Section IV wings of constant-lift distribution (w = constant on 
wing) will be treated. In Section V wings of mixed type are treated. 

These are wings where the pressure is prescribed on some region of the 
wing and the shape on some other region. It will be assumed that the 
second region is flat and at zero angle of attack. In the first, w will 
be assumed to be constant. Only the lifting case is of interest. These 
wings give essentially an interference effect, the interference of a lift- 
ing element on a flat plate. Finally, Section VI treats another interfer- 
ence effect: that of a flat lifting wing on an adjacent flat plate. 

Conical wings will also be classified according to planform. In all 
cases of interest the planform is bounded by two rays from the origin 
which may be considered to be the edges of an infinite triangle. These 
are 

Planform I, in which both edges are outside the Mach cone. 

Planform II, in which one edge is outside, one inside. 

Planform III, in which both edges are inside. 

For each of these cases subcases (a) and (b) are distinguished: 

(a) Both edges are leading. 

(b) One edge is leading, one trailing. 

There is also an intermediate subcase (ab): 

(ab) One edge is leading; the other is an unyawed side edge, 
i.e. parallel to the direction of the free stream. 

For the boundary conditions at the subsonic edges, the statements in 
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Section I-B apply. 

Consider now a supersonic leading edge with a sweepback angle <p . The 
sign of <f> is defined arbitrarily so that <f> is positive if the straight line 
passing through the edge passes through the first quadrant (xj. x 3 0) of 
the plane of the wing. In Figure 2, $ is positive for the right-hand edge 
and negative for the other. In all cases treated in this report the region 



FIGURE 2. SUPERSONIC LEADING EDGES 
of the wing outside the Mach cone from the origin will be flat since here 
w = constant is equivalent to v = constant. From the edge two plane Mach 
waves originate. In the region between these waves and the Mach cone from 
the origin, exact sweepback theory is valid; the velocity components are 
constant and given by 
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. - a Wo,* tan <p 

u — w tan 9 = ■ — 

s/\ - tan^ f 


V = -a. w„ 


(K w, 


w = 


-v 


(1.44a) 

(1.44b) 

(1.44c) 


vl - tan^ <9 Vl - tan 2 9 

It should be remembered that the free stream Mach number is assumed to be 

equal to FT. 

The inclination of the leading edge Mach wave itself is given by 
sweepback theory and is such that its plane is tangent to the Mach cone 
from the leading edge. 

Across the Mach cone a continuous change in the values of u, v, w is 
assumed. Hence on the Mach cone the perturbation velocities are given by 
Equations (1.44) for the part of the cone in contact with the above- 
mentioned regions between Mach cone and Mach waves. On the part of the 
Mach cone in contact with undisturbed fluid u, v, and w = 0. This condi- 
tion holds true in particular for the entire Mach cone if the wing is 
completely inside the Mach cone. The boundary conditions above yield ob- 
vious boundary conditions in the ^ -plane. Consider Figure 3. The region 
between the upper surface, the Mach wave from the right-hand edge, and the 
Mach cone corresponds to the region ACD in Figure 3. The angle Y between 
OD and OA is determined by the formula 

cos Y - tan <(> (1.4 5) 

This formula is easily derived from the fact that the trace of the Mach 
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FIGURE 3. BOUNDARY CONDITIONS IN t -PLANE 
wave in the t, -plane (namely, AD) has to be tangent to the Mach circle. 
Thus if (X is the angle of attack of the upper surface, u, v, and w are 
given by Equations (1.44) on the Mach circle for 0 ±. 0 £ . For 0 > Y , 

the perturbation velocities are zero till 0 is large enough to bring the 
point £ = e under the influence of the left-hand wedge. 

Since 0 is unchanged by the T schaplygin transformation, relation 
(1.45) is also valid in the £ -plane. Subsonic edges are represented as 
points inside the Mach circle on the X-axis, both in the e -plane and in 
the -plane. Condition (1.13 1 ) of Section I-B may be stated when applied, 
say, to w: At a trailing edge (which is represented by a point in the 
£ -plane) the analytic function W has a zero; at a leading edge £ = & , 
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this function may be written W = yv ^ - £ j times a function analytic in 
at 6^ . 

In the 6 -plane, however, one may make an essential reformulation of 
some boundary conditions. In particular the flatness condition (X - con- 
stant may be given a very convenient formulation (Cf. Ref. 4): 


On a flat part of the wing 0w*/#x = 0 


(1.46) 


Consider a direct elementary proof of condition (1.46): Since CL is 
constant, v is constant on the wing and ^v/^x-j = 0. Hence by irrotation- 
ality £w/fl x^ = 0 on the wing in the ^ -plane. Changing over to the €. - 
plane, one observes that on the X-axis, the x^ -direction coincides with 
the y-direction. Thus it follows from the symmetry of the Tschaplygin 
transformation that in the € -plane the curves xj = constant cut the X-axis 
orthogonally. Hence ^w/^x^ = 0 implies ?w/2y = 0. From the Cauchy- 
Riemann equations 9 w*/?x = 0. 

From Equation (1.46) Busemann concludes that w* r constant and that 
this constant then may be put equal to zero since w* is determined only 
within a constant anyway. However, if there are singularities in W, w* 
may take on different constant values on different segments. In particu- 
lar, w* might take on different values on the top and bottom sides of a 
flat wing. This situation may be clarified by considering formula (1.42h), 
which shows that for real er , i.e., in the plane of the wing, the downwash 


is connected with w* by 
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Equation (1.47) furnishes a direct proof of Equation (1.46). Consider 
an edge when e + 0 and a path of integration consisting of a small circle 
with the edge as center. Integrate Equation (1.42b) along this path from 
one surface to the other and let the subscripts u and 1 denote upper and 
lower surfaces. Then 

e 2 -l 

v u - v r +i (w u - > a- 48 *) 

and taking the real part 

, 6 2 -1 . . 

v u " V J = -(' ~ 2 ' fe ) ( w u * - W£ *) (1.48b) 

Thus in the lifting case w* has the same value on the upper and lower 
surfaces. In the symmetrical case the difference (w*^ - w*^) is ZAw^ 

divided by the real constant (1 -e 2 )/2e, where A is the half angle of the 
profile. Equation (1.48b) also shows that when 0, w* u = w£ in order 
to keep the difference v u “ ^ finite or zero. 

Another possible singularity of v and w* is a discontinuity in slope 
on the surface of the wing. Section VI treats of wings where (X takes on 
two constant values on the upper surface. Then dv = dw* = 0 on the wing, 
but at the discontinuity w* has a jump which is obtained from the jump in 
v with the aid of Equation (1.47). 

The discussion just given shows that in the following important case 
Busemann’s statement is unconditionally true: 

On a flat wing of zero thickness and no discontinuity 

in shape, one may put w* = 0 everywhere on the wing. (1.46*) 
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One of Busemann’s fundamental achievements was to point out the appli- 
cation of Equation (1.46 , ). 

L. Oblique Transformation in £ - and 6 -Planes: Equivalence to a Homo- 
graphic Transformation 

In conical flow only the rays from the origin but not the individual 
points are of importance. Each ray (cx^, cx^, cx^) where c runs through 
the positive numbers may be represented by the point 

(±. 5 .. 1 ) 

V X 3 X 3 / 

in the ^ -plane. The oblique transformation induces a certain transforma- 
tion among the rays or among the points in the ^ -plane. Apply Equation 
(1.17) to a point P = (xp x , 1) in the £ -plane. P is transformed into 

w 

the point 


P* 


rx x + a 


1 + ax. 


This point is on a ray which intersects the 2, -plane in the point 

j) 

\ 1 + axj 1 + axj / 

Hence in the -plane the oblique transformation (1.17) induces the 


transformation 


* xi + a 

Xj Xj* = —i 


1 + axi 


(1.49a) 


x 2 -?x 2 * - * 2 


\J 1 - a 2 


1 + ax, 


(1.49b) 
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In the applications of the oblique transformation one usually con- 
siders a transformation which takes a ray of angle fi in the plane of the 
wing into the axis of the Mach cone . Such a transformation is obtained by 
putting a = - tan ft in Equation (1.17) . An arbitrary ray of angle t" in 
the same plane is then taken into a ray 7T* where 


tent.- ta ° * • taa P 

1 - tan tan X 

This is the same formula as Equation (1.49a). 

If the sweepback angle 0 is measured as in Section 
(1.50) should be applied to ^(^ - 0)J. Hence 


(1.50) 


I-K, formula 


tan(f - 0*) = tan <2 - A - tan^ 

1 - tan fi tan(— - 0 ) 
2 


But because tan |(TT/2) - 0 J = cot 0 = (l/tan 0), this formula transforms 
into the same form as Equation (1.50) 


tan 0* 


tan 0 - tan /3 
1 - tan ft tan 0 


(1.50-) 


where 0 is the same as in Figure 2 and 0* is the sweepback angle of the trans- 
formed edge. From Equation (1.50) there follows easily 


/l - 


tan Z * = 


vr: 


tan 2 # \/T"- 


tan 2 x 


1 - tan fi tan T, 


(1.51) 


Now each mapping of the ^-plane into itself induces a mapping in the 
e -plane in the following way: If P is mapped on Fj in the L, -plane , then 

the induced mapping maps Ts(P) on Ts(P 1 ). Since Equation (1.1) is invar- 
iant under the Lorentz transformation, the Laplace Equation (1.31) must be 
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invariant under the mapping which the oblique transformation induces in 
the d- -plane. This induced mapping is then conformal. It leaves unit 
circle and X-axis invariant and preserves direction on the X-axis. This 
invariance follows from the fact that both the oblique transformation and 
the Tschaplygin transformation possess these two properties. Furthermore, 
since in the ^ -plane the origin is mapped on a, in the €r -plane Ts(0) = 0 
is mapped on 


Ts(a) 


A = 


1 - y/l - a 2 

a 


From the theory of analytic functions it then follows that the induced 
mapping must be the homographic transformation 




_ A + e 
= 1 + At 


(1.52) 


The determination of this function is actually a very simple application 
of two principles to be discussed in Sections II, III, and V, namely, ex- 
tension of a function by reflection (in this case the unit circle) and 
determination of its analytical expression from its zeros, poles, and 
Riemann surface. 

The following theorem has thus been proved: 

A Lorentz transformation in the ^ -plane of the type 
(1.49), which takes the origin into the point = & , 
corresponds to a homographic transformation (1.52) in 
the € -plane which takes the origin into 


A = 



a 


(1.53) 
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Expressed in formulas, theorem (1.53) states that if 

^ = Xj + ix£ ( a ) 


‘ f - 1-^*1 * *) * * ' / ^ r « 2 ] 


«= = Ts (£) 

A = Ts(a) 


e * = 


6 + A 


then 


1 + Ae 

e* = Ts( £*) 


(b) 

(c) 

(d) 

(e) 


As a corollary, if Equation (1.52) is applied to a real € ( 6- 
tains the formula 


= T), 


one ob- 


T + A . 
1 + TA' 


1 + at - \/ (l - t 2 )(l 


t + a 


(1.54) 


In the application it is sometimes convenient to use the algebraic 
identity (Cf. Eq. 1.37e): 

_ t + a , 

The right-hand side of Equation (1.54) = j ~ ~ ~ _ a 2) (1.54') 

For the benefit of the reader who distrusts an abstract proof, a 
computational proof of theorem (1.53) is supplied here. It is straight- 
forward but longer than the proof just given. 

Substituting Equation (a) in Equation (b) and multiplying the numer- 
ator and denominator by 

Kl 
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u 

From Equations (c) and (d) 


==- + 2 a 


(l* k) 


| =(j + 6 )> a = 


2 A 


A 2 + 1 


Hence 


.11 - v 2A , — 1.1 > 1 - A 2 , l 1 

(^ + 1+ ^ + 6) + TT7 3(6 +i -)( f + ^) + rrX? ( ¥ + a- --j) 

(<- + -)(e +=-) + i + a2 ( J+J+ £ + e ) 

_ (1 + A 2 ) [e + i + £ ( £ + £ )] + 2A (1 + 6 £) 2 + (1 - A 2 ) [ - <?+ 6g( ■£ -i)] 

(1 + A 2 )(l + £ 6) 1 + 2A [e+ £ + 6. i( £ + £)] 

_ (1 + A z )( £ + 6 ) + 2A(1 + 6 £) + (1 - A 2 )( £ - i ) 

(1 + A 2 )(l + £ £) + 2A( £ + £) 


2 £ + 2A 2 £ + 2A(1 + £ £) 

1 + A 2 + £ i+ A 2 £ i+ 2A( e + e) 


2(1 + & A)( e + A) 


(1 + £ A)(l + e A) + ( e +A)( £ + A) 


2 


£ + A 
1 + e A 


2 6* 


<£ + A _ £ + A 

1 + 1 + fcA 1 + £ A 


1 + £ * e* 


Combining this formula with Equation (1.34) proves that £* * Ts(^*). 

Let (ff and Y be defined as in Section I-K and related by Equation (1,45). 
A Lorentz transformation with parameter a = - tan^ = - b transforms (f> 
into an angle <f> * defined by Equation (1.50 ). The angle Y * s transformed 
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into an angle Y* which, by the previous reasoning, must satisfy the 
relations 


* 


tan <()* = cos^ 

(1.55a) 

iljj* - B 

~ 1 - Be iy 

(1.55b) 


where 


B = Ts(b) 



The following relations may be proved easily and will be of use in 
Section III-G: 


cos 




1 j- cos 
1 - b cosY 


sin 


ll 


1 + b 


' 1 - cos Y 

1 - b cos W 


(1.56a) 


(1.56b) 


In Section I-F methods for generating solutions by oblique transfor- 
mation were discussed. Theorem (1.53) suggests that there should be 
analogous ways of generating solutions in the t- -plane by a homographic 
transformation. 

Although the velocity potential itself does not appear as a function 
analytic in e, there is still an analogue of the method of generating so- 
lutions by transformation of the velocity potential. Let U, V, and W be 
complex functions constituting a solution of some conical flow problem. 

In particular they have to satisfy the compatibility relations (1.42). 

Then the following three functions also satisfy relations (1.42) and 
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constitute a solution for another conical problem: 

U a (e) = ■ ) + A * U(01) + ^ 5 - W(7?) (1.57a) 

1 “ A* 1 - Ar 

V*(e) = V(7l) (1.57b) 

W*(fe) = 2A U(^)+~ - + ~ A y W (71) (1.57c) 

1 - A 2 1 1 “ A^ 

where 

77 = ±±A 

L 1 + At 

This conclusion follows from Equations (1.21c) and (1.3 7) and may also 
be checked by direct computation. 

Using the method expressed in Equations (1.22), one may also put 

w a (e) = W ( 77 ) (1.58) 

It will be seen from examples to be given later (Cf., e.g., planform 
lib in Section III-H) that a general formula like Equation (1.57) for con- 
structing U a and V" is not possible. Hence these functions will have to 
be computed from with the aid of Equation (1.42) in each separate case. 

In the applications of Equation (1.58) the following question often 
arises: If W is the solution for a flat wing, is the same true for W a ? 

The hypothesis implies that 

dv = dw* = 0 on the wing (a) 

From condition (a) and Equation (1.58) there follows: 

dw a * = dv = 0 on the corresponding wing (b) 

Thus the second wing is flat (v = constant) unless V has a singularity 
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there. If W is regular on the first wing, W a is regular on the second 
wing. Then the only possible singularity of V is at e = 0. The nature of 
this singularity is given by Equation (1.42c). Thus the second wing is 
flat if it does not contain the origin. 

M. Prandtl-Glauert Transformation Reduction to Mach Number 

In the present report it is generally assumed that the free stream 
Mach number M is equal to "fz; It will now be shown that within the line- 
arized theory a wing problem at any other supersonic free stream Mach num- 
ber may be reduced to an equivalent problem at M = \flT. 

Let the original problem be to find the perturbation velocity caused 
by a wing placed in a flow of free stream Mach number M. Denote by S m the 
three-dimensional space in which this flow takes place. The problem will 
be reduced to the study of an equivalent wing placed in a free stream Mach 
number = > T~2 in a three-dimensional space S. (S and S m are two copies of 
ordinary infinite three-dimensional Euclidean space.) A certain corre- 
spondence or mapping between the points in S and S m will be set up which 
will in particular define the equivalent wing: It consists of the points 
in S on which the points of the wing in S m are mapped. From the properties 
of this mapping there will follow a procedure for constructing the solution 
for the flow in S m from that for the flow in S. 

The mapping is the well-known Prandtl-Glauert transformation. A 
point P in S™ will be mapped on the corresponding point P* in S by the 
rule 

P = (xj, x 2 , x 3 ) >P* = (x lt x 2 , mx 3 ) (1.59) 
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where 

1 

m = tan M. = - — 

/m 2 - 1 

This correspondence between points will induce a correspondence be- 
tween functions. Two functions, f in S and in S m , are corresponding 
functions if 

^(P) = f(P*) ( 1 . 60 ) 

Similarly, the two differential operators 0 and 0 m will correspond to each 
other if 


value of operator 0 m applied to f m at P equals value of 
operator 0 applied to f at P*. (1.61) 


Consider an example: 

If 0 is —— , then 0™ = - y|— (1.62) 

This correspondence may be proved as follows: 

^(XJ. x 2 , x 3 Hx) - f m (x lP x 2 , x 3 ) 

lim 

Ax-*0 Ax 


-(— ) = L 
m\3x-/ m 


- lim 
m Ax 0 


f(x^x£, mx^+mAx) X£, 1 * 1 x 3 ) 

m A x 



The two wave operators Q and Q m in the two spaces are defined by 


9xj 2 8 x 3 2 3x3^ 
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2 2 

Q m = — + — — 

0x^ 2 9 x 2 2 




2 


By the method used immediately above it is proved that 

Q corresponds to Q 111 (1.63) 

The definition (1.61) and the theorem (1.63) imply that: 

If Qf = 0 at all points in S, then C^f 111 = 0 at all 

points in S 111 . (1. 64 ) 

In other words, if f is a solution of Equation (1.1) , then f m is a 
solution of Equation (1.1*). In this way solutions to Equation (l.l 1 ) may 
be generated from solution (1.1) by the rules (1.59) and (1.60). This 
procedure is very similar to the previously studied generation of function 
by the transformations mentioned in Section I-C. In the present case, 
however. Equation (1.1) is not invariant under the transformation consid- 
ered, and for this reason the generated function f 111 will not satisfy 
Equation (1.1) but the transformed Equation (l.l 1 ) instead. 

In order to make use of relation (1.64), one has to investigate what 
boundary conditions the transformed function satisfies. Apply Equation 
(1.60) to the velocity potential <f) , forming a potential for the flow in 
S™. Form the velocity components u, u m , v, etc. from these two potentials. 
One obtains the following relations (Cf. Eq. 1.62): 
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* m (P)= f(P*) 
u m (P) = u(P*) 
v m (P) = v(P*) 
w m (P) = mw(P*) 


JPL 


(1.65) 


Now let there be given a wing in S m of prescribed shape, satisfying 
Equation (1.6). If Equation (1.59) is applied literally to the wing, the 
equivalent wing has a different angle of attack because of the affine dis- 
tortion. However, as explained in Section I-B, in the mathematical treat- 
ment the wing is replaced by its projection on the XjX^-space, which is 

given a certain distribution of OL , i.e., of v. It is to this projected 

i 

region that Equation (1.59) should be applied in forming the equivalent 
wing. It then follows from Equations (1.64) and (1.65) that Equation 
(1.65) is the correct rule for generating a solution for the original wing 
from that of the equivalent, provided the equivalent wing is prescribed to 
have the same angle of attack as the original wing at corresponding points. 

An alternate scheme is used if the boundary conditions are given on w 
instead of on v. In this case let w have the same value at corresponding 
points of the two wings. Then one obtains the functions etc. for 

the equivalent wing, forms the functions etc. as in Equation (1.65), 

and then multiplies them by the quantity m. 

I 

If one has mixed boundary conditions (Cf. Section I-B), the first of 
the above shcemes works if the only value prescribed for w is zero. Simi- 
larly, the second scheme works if the only value prescribed for v is zero. 
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II. FLAT SYMMETRICAL WINGS 
This section discusses flat conical wings whose profile (i.e., inter- 
section of wing and the plane = constant) is a wedge of half angle A 
(Cf. Fig. 4). The wings are assumed to be symmetrical relative to the 



FIGURE 4. BASIC CASES OF FLAT SYMMETRICAL 
WINGS 

plane of the wing. The general characteristics of symmetrical wings were 
discussed in Section I-B. In the plane of the wing the boundary condi- 
tions are 

v = Aw a on the upper surface of the wing 

v = - A w^ on the lower surface of the wing (2.1) 

v = 0 off the wing 
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On the Mach cone the values of v are given as described in Section I-K. 

As for boundary conditions at an edge, consider a subsonic edge, 
leading or trailing. In the e -plane such an edge is represented by a 
point (-/ <6 1 <l ) and the boundary conditions in the neighborhood are as 
in Figure 5. 


Y-axis 


Wing 


A v=v o= 


Aw e, 


=Edge of Wing 

r 


, v 0 =-Aw 00 


v = 0 


X -axis 


FIGURE 5. SUBSONIC EDGE OF SYMMETRICAL WING 
It is seen that the given values of v are v o /tT times the polar angle 0 
with apex at £ ^ ; namely, if 0 = 0, + TT , or - 7T , then v = 0, v Q , and -v Q , 
respectively. The analytical expression for the fact that v is propor- 
tional to 0 is that it is the imaginary part of a logarithmic singularity 
at . More precisely 

✓ 

v 

V=-i^p-in(£ - 6 X ) + i(e) (2.2) 


where f(e) is analytic at and purely imaginary along the X-axis. Thus 
f (t) does not contribute anything to v = Re(V) there. 

From Equation (1.42b) it follows that 


dW V o 


d€ IT £ -1 


2 ^ 


V *(* -!) 


f’(e) 


(2.3) 
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Hence for 6 t 0, W also has a logarithmic singularity. However, w is the 
real part of the logarithmic singularity and hence becomes logarithmically 
infinite at & - €r j , whereas w* behaves qualitatively like v. In the ex- 
ceptional case where the side edge is along the axis of the Mach cone, 

/ =0 and V still has the same type of singularity but W is regular. 

Although the function (2.2) fits the boundary conditions in Figure 5, 
there are other possible solutions. These will have to be ruled out by 
special considerations. This subject will be discussed in Section II-B. 

As pointed out in Section I-B, the fact that there are sufficient 
boundary conditions for one perturbation velocity makes the problem of 
finding solutions for symmetrical planar wings very easy. Another circum- 
stance which simplifies the problem is that there is only pressure inter- 
ference but no downwash interference between two such wings. Hence the 
boundary conditions are not affected by the interference. Thus in the 
present case it is sufficient to obtain the solution for planforms Ilab 
and Illab (Cf. Fig. 4). The other solutions for planforms I, II, and III 
may then be obtained by superposition. In particular the solution dis- 
cussed in Reference 17 may be derived in this way. 

A. Planform Ilab 

The boundary conditions are shown in Figure 6. As has just been ex- 
plained, the singularity at the edge may be expressed by the function 


■1 v r 


rr 


Jne 


( a ) 


This function gives the correct values of v on the X-axis provided the 
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JPT 


fcr -plane is cut along the negative real X-axis. The two discontinuities on 


the Mach Circle may also be obtained as logarithmic singularities. If e 



FIGURE 6. BOUNDARY CONDITION IN e -PLANE 

varies in the positive direction on a small half circle inside the Mach 
circle (Cf. Fig. 6), then v decreases by v Q . This type of discontinuity 
suggests that V is of the form 

iv o 

In (e- e^) + function analytic at ^ (b) 

Combining functions (a), (b) , and the corresponding function for the dis- 
continuity at 6 ^ gives 

v = [-Me - e x )(e - (c) 

Actually the function defined by (c) satisfies all the boundary conditions. 
Since 

Re [-iXn(€- fr )] =arg(£-e 1 ) 
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and on the real axis arg (e - e j) = - arg(fr - ej , the first term gives 
no contribution to v on the real axis , and the second term gives the 
correct boundary conditions. That conditions on the Mach circle are 
also satisfied will follow later on (Cf. proof of Eq. 2.8) when v is 
evaluated in physical coordinates (r , 6) from Equation (c). 

But first U and W will be evaluated from Equations (c) and (1.42). 
From Equation (c) 


dV 
de “ 



6 -6 



(<*) 


From relation (d) and Equation (1.42b) 

dW 2v 0 e 

de ir e 2_i 


6 




Since & 

1 




(e) 

1 


^ ^ = ^isin ^ 


(f) 


Combining relations (e) and (f) , 


IV 


W = 


ir siny 


Jt n 


6 - 6 - 

t 

6T7" 


+ c 


(g) 


To evaluate the constant of integration c, one observes that. 


according to (1.44) and (1.45), 

~ v o 

at 6 = -1, Re(W) = w = w = — : — - 

v ' o sinf^ 
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But as may be seen by elementary geometry 


Hence 


Re 


7T W 

w _ 


-i in ~ 1 


-1 - e- 


e - £ 


= -iin 


£ - e 




+¥+ -rr 


(h) 


Similarly U may be evaluated 


dU 1 , 1, dW -i v 0 Iv 

de“ ’2 ( e+ ^de ~ 2 7T sin'f 


-1 


6 - C 'l e - e 2 


(i) 


After some algebraic manipulations this formula becomes 


dU 
d e 


"VV x: 

7T 


COS 


y. 1 


cos 


v _ . 


sin y €• - sin €■-€-. 


e 


(j) 


Integrating formula (j), using Equations (1.44) and (1.45), and determining 
the constant as for relation (g) 


U = 


7T 


- iin TT" g 2 - +y+T 

1 


“ ^ inf 

qr 


(k) 


Next u, v, and w will be evaluated as functions of the physical coordi- 

iQ y 

nates r and 0. Consider an arbitrary point re in the £ -plane and the 

id 

corresponding point Re in the £ -plane. 

Using the notation in Figure 7 


Re( - iin(<rj - C-) = arg^ - e) - 0 1 
Re( - iin(e 2 - e) = arg(e 2 -e) = 0 2 
These angles are connected with R and 0 by 


(/> 
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tan 


sin 


f - R 


sin@ 


1 cos y - R cosG 


tan <p ? = 


sin - R sin0 
cos - R cose 



(m) 


FIGURE 7. NOTATION FOR EVALUATION OF u, v, AND w IN € -PLANE 
From relations (m) 


, /a t ± x Z R sine (R cos e - cosff 1 ') , . 

(^i + ^ 2 ) = 1 + RZ - 2 RZ s in z 0 - 2 R cosG cos Tp 


tan( ! + f z ) = L ( 1 cos9 - r cos ^ ] (n>) 

^ " (1 “ >/l " r^) sin^e - r cos e cos ^ 


In a similar way one may derive 


cos (♦2 -*i +f+r) 


r cos 6 - cos 


\/(^ r cose cos y y 2 - r^ sin^e sin^' 
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The results derived are summarized in Equations (2.4) to (2.7): 

For a wing of planform Ilab, sweepback 0 , and symmetrical wedge profile 
of half angle A, the perturbation velocities are as follows. On the wing 


outside the Mach cone 


u = u Q = - w Q tan </> = 


A Wqq tan (j> 
«/ 1 - tan^ <P 


v Q = A on the upper surface 


v = -v Q = - A w^ on the lower surface 


w = w_ = 
o 


~ v o . -Aw^ 

1 - tan^ f fl - tan 2 <p 


(2.4a) 


(2.4b) 


(2.4c) 


Inside the Mach cone the complex expressions for the velocity are 


U= f-iin 6 

T I * J IT 

= ^?-|[in(6 - 6-^(6 - e 2 ) n 6-J 


*■ IF - U “* 


(2.5a) 


(2.5b) 


(2.5c) 


Here cos Y = tan <£ , e ^ = e , 6^ = 

Putting arg(^ -6) = ar g(^2 ~ e ) = $ Z (Ct m Fig. 7), and 

i0 

6 = R * e , the real parts of the functions above may be written 


"= £[v«i + v ' +r ] • £ in 


(2.6a) 
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v - V ° 

v ~ ir 

® - ( ^ 1 + 

( 2 . 6 b) 


■Wo 

W = j 

IT | 

^» 2 + tr] 

( 2 . 6 c) 


With the aid of relations (n 1 ) and (o) the values of the brackets in 
Equations (2.6) may be given as a function of the physical coordinates r 
and 0 


cos(<f> 2 -p l + +TT) 


r cos 9 - cos ^ 


J (1 - r cos 8 cos ^ - r 2 sin 2 6 sin^ 


(2.7a) 


tan( 4 >j + <t> 2 ) = 


sin 8 [ (1 - \/l - r^) cos 8 — r cos J 
1 - a - r^) sin^0 - r cos 0 cosV 7 


(2.7b) 


For points in the plane of the wing Equations (2.6) and (2.7) may be 
simplified considerably. Use is thereby made of the fact that 0 is either 
0,7 T, or -7T and hence sin0 = 0 and r cos 0 = t (Cf. Eq. 1.36). 

In the plane of the wing 


u ° t - tan $ 

u - — arc cos — ~ 

7 T 1 - t-tan <p 


v o f 1 -/ 

ir /n 1 1 1 


l 


V = _ e 

7 r 




t - tan <f) 


w = — arc cos 1 . t . tan 0 


(2.8a) 

(2.8b) 

(2.8c) 


Finally on the Mach cone Equation (2.7b) reduces to tan 0 = tan (<P f + • 

Hence 0 - = n 7T where n is an integer. At the point £ = 1 , n 

must be zero by Equation (2.8b). Hence by analytical continuation it re* 


mains zero on the Mach circle up to singularities and C- . Here the 

J 
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jump is found by the method described at the beginning of Section II. In 
this way it is seen that Equation (2.5b) actually satisfies the boundary 
conditions on the Mach circle. This fact could also have been seen geo- 
metrically from Equation (2.6b). 

B. Planform Mab. Principle of Reflection 

The boundary conditions as illustrated in Figure 8 are very similar to 



FIGURE 8. SYMMETRICAL WING: PLANFORM IHab 
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those in Section II. The main difference is that the logarithmic singu- 
larities are all on the real axis of the e -plane. Since the wing is en- 
tirely within the Mach cone, there are no singularities on the Mach circle; 
here u = v = w = 0. The function V has two logarithmic singularities in- 
side the Mach circle, at € = 0 and ^ = B, which points correspond to the 
two edges. However, these two singularities are not sufficient to deter- 
mine V. Some information about the behavior of V in the rest of the 
€: -plane will be necessary and can be secured with the aid of the principle 
of reflection. 

A proof and exact formulation of this principle may be found in text- 
books on analytic functions, especially those with the geometric (Riemann) 
approach (Cf. p. 218 of Ref. 18 or p. 372 of Ref. 19). The principle may 
be conveniently formulated in geometrical terms. One may consider the 
solution V (e) as a mapping of the e -plane into the V-plane; to each point 
x + iy in the 6 -plane there corresponds a point V(e) = v + iv* in the V- 
plane. In this way a straight line or a curve in the e-plane will corre- 
spond to some curve or straight line in the V-plane. 

The principle of reflection (Cf. Fig. 9) indicates that if a line 
is 'mapped analytically on a line and the mapping is defined on one side 
of L^, then it may be extended to the other side of Lj in such a way that 
points which are symmetrically located (images) with respect to L,\ are 
mapped on points which are symmetrically located with respect to L^. 

In this statement any one of the straight lines may be replaced by a 
circle if by image of a point with respect to a circle is meant the recip- 
rocal point. Two points are reciprocal with respect to a circle if they 
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FIGURE 9. PRINCIPLE OF REFLECTION 


are on the same ray from the center of the circle and the product of their 
distances to this center is equal to the square of the radius. 

The principle of reflection will be of fundamental importance for 
many solutions given in this report. It will now be applied to the present 
example. On the unit circle in the C -plane, v = 0; in other words, V is 
purely imaginary. Expressed in geometrical terms this fact means that the 
unit circle in the £ -plane is mapped on the imaginary axis in the V- plane. 
That is the image of ^2 * n *h e un ^ circle means algebraically that 

• Similarly, if is the image of V 2 with respect to the imagi- 
nary axis, then Vj = -V^. Hence V must satisfy the relation 


V(l/e) = -V(£-) (2.9) 

More precisely it is known that V is analytic inside the unit circle ex- 
cept for the slit OE (Cf. Fig. 8). Since v = 0 on the circumference of 
the unit circle, V may be extended to the outside of the unit circle by 
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the rule (2.9), and it will be analytic there except at the slit which is 
the image of OE. The reciprocal of the point E: €- = B is the point 
E*: € = l/B. The reciprocal of the origin is the point at infinity. Hence 
the image of OE is the slit extending along the real axis from l/B to in- 
finity. From Equation (2.9) there follows that v = - on top of the slit 
and v = + v Q on the lower side. Thus there is the same type of logarithmic 
singularity as at the edge E. One may think of E* as the edge of the re- 
flected wing. When €* varies counterclockwise on a small circle around E*, 
V increases by 2v Q . This fact determines the strength of the singularity. 
Hence all singularities of V are known, and the solution may be written 
down immediately as for planform Ilab: 


V = 


"i vq 
7T 



e - B)(B6 



v 


o 


= A w 


00 


( 2 . 10 ) 


where \ = half angle of profile. The functions U and W are then easily 
computed from Equations (2.10) and (1.42): 




(2.11a) 


o 

tF 


(2.11b) 

That the constant of integration has been chosen correctly may be 
seen from the fact that if & = 1, Re(U) = Re(W) = 0 which is the correct 
boundary v* ue. It is interesting to note as an additional check that 
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- (e - B)/(B6- -1) is the homographic transformation previously encountered 

in a different context (Cf. Section I L) . This transformation maps the unit 
circle of the 6 -plane into the unit circle of the -plane. Hence when e 
is on the unit circle, the arguments of the logarithmic functions in Equa- 
tion (2.11) are applied to numbers of absolute value = unity, and u and w 
vanish on the entire unit circle. 

For an arbitrary point re 1 ^ in the 4 -plane, w has the value 

- B 


v rt 2B q_ 

W = ^P rr&r STl 


Re 1 
BReiS - 1 


(2.12a) 


where (Cf. Eqs. 1.29 and 1.37h) 

2B b 

v o= Aw oo> R = Ts(r), B=Ts(b), and i - ~ j ~ _ b2 

i0 . 

In the plane of the wing 6 = Re = T = Ts(t). Hence 


Vo 2B 

0 

T - B 

v o b 

0 

1-bt- \/ (1— b 2 )(l — t 2 ) 

W — __ 

IT i - B 2 

X r 

BT - 1 

b^ 

-X n 

t - b 


(2.13a) 


Here use has been made of Equation (1.54). Similarly 

>2 


u = 


Vo 

TT 


( , - R -r 


B 2 




Re 10 - B 


BRe ie - 1 


) 


(2.12b) 


and in the plane of the wing 


u = V o 


Jn l T l - 7= f 

TT \ V 1 - b 2 


■it 


1 - bt - /(I - b Z )(l — t 2 ) 


t - b 


(2.13b) 
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C. Solutions with Higher-Order Singularities: Uniqueness 

Some comments will now be made about the uniqueness of the solutions 
obtained in Sections II-A and B. 

Consider, e.g., planform Ilab. Denote by Vj the solution obtained in 
Section II -B and by another solution with the same boundary condition. 
The difference = Vj - then satisfies the following boundary condition: 


v 3 = Re(V 3 ) 


0 on real axis and on unit circle 


(2.14) 


If the function V 3 is to be regular in the interior of the unit cir- 
cle, condition (2.14) implies that it is identically zero (or equals an 
imaginary constant). In this case Vj = except for irrelevant imaginary 
constants. 

Nonvanishing solutions of condition (2.14) may of course be had by 
assuming various unnatural singularities of V^. However, assume that the 
only singularity is at the origin £ = 0 and that except for this point 
is regular in the unit circle. 

Then by the principle of reflection V 3 is regular on the whole Riemann 
sphere ( e -plane and point at infinity) except for the singularity at e = 0 
and the reflected singularity at ^ = 00 . If the singularity is a pole of 
order n, the function must then have the form 

v 3 = V< n > (e) = ki(t n +^ 1 ) (2.15) 

where k is a real constant. Hence 


v 3 = He(V3> 


= k^n - R nj 


(2.16) 


sin n 0 
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The corresponding W -functions are given by 

(e) = 2k in& (2.17a) 

W< 2) (*) = 4k (e - \) (2.17b) 

(n) . 

Thus for n > 1 the pressure distribution corresponding to V gives 
an infinite force on a finite area, and the solution may be rejected be- 
cause of this fact. However, the force corresponding to the solution 
is finite since only has a logarithmic singularity at the edge. Thus 

the only way to reject this solution is to make explicit use of the bound- 
ary condition (1.13 "). 

If one allows discontinuities across slits in the unit circle, poles 
of fractional order may be used. Let the wing correspond to a slit in the 
e -plane. Consider the function 


w = i(-(?+ 


= (_i — - ’/r) sin 

Vr 


(2.18) 

(2.19) 


The corresponding downwash function is 

(2 . 20) 

6 


v = - 7=) cos | - ^ (/r 3 -^ cos f < 2 - 21 ) 

Evidently v vanishes on the wing 0 = -f 7T and on the unit circle 
r = l/R. The pressure is infinite at the leading edge but gives a finite 
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force for finite areas. Thus again it is possible to add this solution to 
the solution of Section II-B unless boundary condition (1.13 11 ) is used. 

It is interesting to note that Equations (2.18) and (2.20) have the 
symmetry properties of the lifting case. Thus it may be seen that in 
spite of the fact that the equation and the boundary conditions are anti- 
symmetric in v, the resulting solution does not necessarily have this 
property unless uniqueness may be proved (Cf. Section I-D). 

III. FLAT LIFTING WINGS OF ZERO THICKNESS 

The general characteristics of flow fields induced by lifting wings 
of zero thickness have been discussed in Section I-B. When the shape of 
the wing is prescribed, finding the solution is difficult because the 
boundary conditions are mixed (involving both v and w). An unpleasant 
characteristic connected with a lifting wing is that it generally induces 
an upwash (positive or negative) in its own plane off the wing. Thus it 
is in general not possible to generate solutions by the very simple method 
that was used in Section II. 

The prescribed shape gives the condition on v. When the surface is 
flat, this condition is v = constant. In this case, which is the only one 
considered in Section III, the condition on v may be eliminated. This im- 
portant method (drawn from Busemann) will be discussed in connection with 
planform Ilab. More general shapes (cambered profiles), both conical and 
nonconical, may be constructed from the wings introduced in Section VI. 

Especially simple are those solutions which may be obtained from the 
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solutions in Section II with the aid of the general principle to be described 
next. 

A. Cases Reducible to Symmetrical Flow 

Consider a conical symmetrical wing which may have any camber, but 
whose two edges are leading and supersonic. Then the Mach cone from the 
apex is disconnected by the wing into two separate regions without any 
contact. There is no contact between upper and lower surfaces around the 
edges since these edges are supersonic. Thus the flow above the wing is 
completely independent of the flow below the wing. Hence if the lower 
surface is kept unchanged but the upper surface is modified, the flow be- 
low the wing is unchanged. Modify the upper surface so that the wing 
changes from symmetrical to a wing of zero thickness. Then below the wing 
the flow is known from the symmetrical case. Next one obtains the flow 
above by the general symmetry rules for the lifting case given in Section 
I-B. 

B. Planform la 

This is the wide delta wing. Both edges are outside the Mach cone 
and leading. Hence the solution is obtained from Section II by the general 
principle just discussed. Let the wing have an angle of attack = (X . Con- 
sider the solution for the symmetrical case where the half angle of the 
profile is = A . Then the perturbation velocities below the wing in the 
lifting case are obtained from the formulas for the symmetrical case by 
replacing A by (X . To obtain the flow above the wing, A will have to be 
replaced by - (X . 
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C . Planform lb 

Here one edge is trailing and the wing is completely outside the Mach 
cone . In the applications such a case may occur at a wing tip if the rake 
angle is larger than the Mach angle. This case was discussed briefly by 
Busemann in Reference 4. Conditions on the wing are entirely supersonic 
and obtained by sweepback theory. For this special planform a flat wing 
is also a wing of constant-lift distribution. Within the Mach cone there 
is no lifting surface but a flow field (vortex field) induced by the lift- 
ing wing outside the cone. The solution is easily obtained by a superpo- 
sition of two wings of constant- lift distribution with planform Ilab (Cf . 

Ref. 4). Such wings will be considered in Section IV, where the superpo- 
sition method will also be discussed. 

D. Planform Ilab: Pressure 

This planform is of importance in considering the tip effect of an 
unyawed rectangular wing and applies also when the sweepback of the leading 
edge is different from zero but smaller than the complement of the Mach 
angle . The planform of the wing and the boundary conditions are illustrated 
in Figure 10. Busemann indicates how this case can be solved by using 
the principle of reflection. (He actually treats the slightly less general 
case <J> = 0). The solution will here be given in some detail in order to 
illustrate the use of certain very important principles. 

As explained in Section I-K , w* may be assumed to be zero on the 
flat wing AO . 

Next use the principle of reflection as explained in Section II-B . 
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FIGURE 10. LIFTING WING: PLANFORM Ilab 
The function W is analytic inside the unit circle although not necessarily 
across the slit AO. The fact that w* = 0 on this slit means in geometri- 
cal language that AO in the £ -plane is mapped on the real axis in the W- 
plane . Hence W may be continued analytically across the slit in such a 
way that points which are images with respect to AO are mapped on points 
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which are images with respectto the line w* = 0 in the W -plane. This 
fact may be expressed algebraically as follows (Cf. formula 2.9): 

w(e) =w(F) (3.1) 

But Equation (3.1) implies that w is symmetric with respect to the 
X-axis, whereas it is known that in the physical problem w is antisymmetric. 
This apparent contradiction is resolved by introducing the Riemann surface. 
Start with a function W analytic in the unit circle except across OA, 
which represents the physical solution and is antisymmetric in W. If this 
solution is now continued across OA by reflection, one arrives at the 
lower sheet of the Riemann surface of the function W. It should be empha- 
sized that any point may be reflected in OA, not only points directly 

% 

above it. Thus the whole unit circle may be reflected in OA, and the 
image will form a second sheet of the Riemann surface which is a unit 
circle like the first sheet. The two sheets are joined along the slit OA. 

In Figure 11 a closed path around the origin is drawn. The function is 
analytic along the entire path. The fact that the path actually closes 
after two laps results from the assumption that w* = 0 both on top and 
below the slit OA. If its valie on the lower side of the slit were a con- 
stant different from zero, reflection would have given an infinite number 
of Riemann surfaces instead of just two. 

The main advantage of extending the function is that the discontinuity 
OA has disappeared. On the Riemann surface the function is analytic across 
OA. Thus there are no boundary conditions in w* or v, but only in w. The 
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Dotted Part of 
Contour is on 
the Lower Sheet 



FIGURE 11. CLOSED PATH AROUND THE ORIGIN ON THE RIEMANN SURFACE 

boundary is now a connected closed curve which represents the circumference 

of the two-unit circle. It runs like the closed curve in Figure 11. The 

values on the part of it which is on the upper sheet are given from physical 

considerations (Cf. Fig. 10). The values on the lower branch are given by 

reflection; e .g. , at the point P in the physical plane w = -w Q but on the 

corresponding point on the branch beneath it w = w^ since that point is 

the image of a point Q in the physical plane where w was prescribed to be 

w . 
o 

This is the first example of Busemann's method for treating mixed 
problems, mentioned in Sections I-B and I-K. In the original mixed problem 
the boundary conditions are given in terms of both w and v. With the aid 
of Equation (1 .45) the condition on v was transformed into a condition w*. 
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Then with the aid of the principle of analytic extension by reflection, 
the condition on w* was replaced by conditions on w on the extended part 
of the Riemann surface. Although the domain obtained after the analytic 
extension is a double -sheeted circular disk, it is actually easier to 
handle than the original disk since it contains no cut across which w is 
discontinuous. 

Now there is a standard way of unwinding such a double -sheeted sur- 
face, i.e., with the aid of the transformation 


V =fe 


(3.2) 


This transformation maps the double -sheeted surface obtained by extending 
the £ -plane into a single-sheeted surface in the V -plane, and the boundary 
value problem can now be formulated in the V -plane. This fact is illus- 
trated in Figure 12. Here is the transform of the point C on the upper 
sheet, the transform of the corresponding point on the lower Surface, 
etc. In particular the wing OA appears in the imaginary axis in the V- 
plane; the segment OAj represents the upper surface and OA^ the lower 
surface. 

W is analytic in the unit circle but has certain singularities on the 
circumference. These are of the same nature as those encountered for plan- 
form Ilab in Section II. Thus the jump in w at the four points V ^ may be 
taken care of by appropriate logarithmic functions. This procedure leads 
to the following solution: 


W(fr) = + 


7T (v-^ 2 )(y_y 3) 


(3.3) 
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V V? as in Figure 12 

w = w o on the top side of the wing outside Mach cone 
+ a 

\xr — j rl L 

o / 7 "~ 

Vl - tan 

It will follow from the evaluation of Equation (3.3) in physical coor- 
dinates that this solution actually satisfies all the boundary conditions. 

It is especially easy to evaluate w on the wing. However, in order to il- 
lustrate the technique, the general expression for w at any point inside 
the Mach cone will be derived. 

The following abbreviation will be used: 
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N = \/~R 

From Figure 13 and formula (3.3): 




(ir 


V -Plane 


All angles 
0. positive 
' if 
- e ~~E 


(3.4) 


~ 4 z) - (^+^3) 

1 -tf > 4 + K • 27tJ 




‘A\ 

-• / L \ 


1 H 


AT 

*^3^ — 



(3.5) 


FIGURE 13. NOTATION USED IN EQUATION (3.5) 


By considering a point on the real axis of the P -plane, it is seen 
that the integer k has to be given the value k = 1, because such a point 
corresponds in the physical plane to a point in the plane of the wing but 
not on the wing. Here w has to be zero because of antisymmetry and con- 
tinuity. Thus 


Tw 


= - <®2 - 0 1 > 


(3.6) 


where 0^ and are explained in Figure 13. Figure 13 also gives the 
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relations 


tan0^ 


sin y/2 - N sin 0/2 
cos Y/Z - N cos 0/2 


, sin y/2 - N sin 0/2 

* an ^2 “ cos y/ 2 + N cos ■072 


^ sin Vy2 + N sin 0/2 

tan ^3 = cost/ 2 + N cos 0/2 


> 


* ^ _ sin V 72 + N sin 0/2 

^4 " cos r/2 - N cos 0/2 

y 


From relation (a): 


- 2 N sin 0/2 (cos */ r /2 — N cos 0/2) 
tan(^ - 0 4 ) = l + N 2 cos 0 - 2 N cos Y~/~2 cos 0/ 2 




/iK j . -2N sin 0/2 (cos'4 r /2 + N cos 0/2) 

tan(0 -0 ) = 1 

3 1 + N 2 cos 0 + 2 N cos t/2 cos 0/2 


y 


From relation (b): 


(a) 


(b) 


tan (p x + p 2 - 0 3 - 0 4 ) = 

2N^cos 0/2cos^/2 sinQ-2N sin0/2cos V / /2-2N^sinQ/2 cos W2 cos 0 

-(4N 2 sin 2 0/2cos 2 V r / 2 )+N 4 sin 2 0+l+N 4 cos 2 0+2N 2 cos 0-4N 2 cos 2 (0/2)cos 2 f/2 


2N sinQ/2 cosV72 (N 2 -!) 

= 2 2N 2 [cos 0 - cos y ] + N 4 -2N 2 +1 

_ 4 sin 9/2 cos / 2 

n r I ^ n 2 -i 

n 2 -i l j n 


(c) 


By Equations (3.4) and (1.37i) 
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(d) 


Also cos y = tan 

cos f/ 2 = y?-0?2(y/2) tl _ J tang + 1 


Hence relation (c) may be written 

. T / , .iv /x . 1 vl sin 58 s/tan 0 + 1 • 277 

[ (♦, * h) - (*s +MJ y_sf (cose . t ,„ 0)t y^L 

Simplifying Equation (f) and comparing with Equation (3.5) give the final 
formula: 


w = 


A woo 




.arc tan * sl “ t Mta rr» + l)r(l - r~ 
1 - r(l - cos 0 + tan <p ) 


(3.7) 


tan <p 

This formula gives the value of w at a point (r, 0) inside the Mach 

cone if the value of w on the upper surface of the wing outside the Mach 

cone is = + w = — j - ^ w 

° \Jl - tan <p 


By elementary trigonometry Equation (3.7) may be transformed to 


w = 


(X w< 


— TT-^-r- arc cos 1 - r(l- cose + tan 4>) _ (3.7*) 

ITv/l tan <P + r^(tan^ <f - sin^0)-2r cos 0 tan (p 


The advantage of the second form (3. 7*) is that on the wing itself 
(0 = IT) it simplifies to 


^w ( 


w = 


TT V 1 " tar * ^ 


arc cos 


1 - r (2 + tan<j>) 
1 + r tan <p 


(3.8) 
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or using the notation of Equations (1.3 5) and (1.36a) 

w = gWp. arc cos 1 (3.8') 

IT f\ - tan 2 1-tanTtan^ 

Equation (3.8*) is easily derived from Equation (3.6) by the above method 
if the cosine function is used instead of the tangent function. 

It should be remembered that this formula is valid on the wing only 
and that hence X takes only negative values. For points in the plane of 
the wing but not on the wing. X is positive and w of course is 0 by anti- 
symmetry. Note that if the range of arc cos is restricted to angles be- 
tween -IT and IT , it is still a two-value function which for the same argument 
gives a positive and a negative value. The positive value corresponds to 
a point on the top surface of the wing. 

For <f> = 0, Equation (3.8*) specializes to the formula given by Buse- 

mann (Cf. Ref. 5). 

E. Planform Ilab: Downwash and Sidewash 

The computations proceed along the same lines as those in Section II. 
First introduce the variable V - V& into relations (1.42). 


dU _ 
dy = 

_ i_ 
" 2 


dW 

dxJ 

( a ) 

dV 

1 1 


dW 

(b) 

dll ~ 

2 1 

) 

d ^ 



Then introduce the functions F and G by 

dF = - _ H . _L_ dW ( c ) 

iw Q JA 2 
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-it 2 

dG =1 — » dW 
o 


The following expansions into partial fractions will be needed: 


l _ j ±_ _L_ . L 

0 .2 .2/ 


-a) a^ a>^ ^ a^(^ -a) 


2 a 2 

y > a 

— = y + a + — — 

V - a V -a 


From Equation (3.3) 


1W 


dW 

dV TT y 


♦^Z* ^ 


where 


and 


+ is used for j = 1,4 


- is used for j = 2,3 

Combining relations (c) , (e), (f) , and (g) 


41 =z. 

dv — 


1 


1 


z (>s-y.) J 


where + and - are used as in relation (g) . 
Since 


i f -i£ 

. 2 2 

= e > ^2 = - e 


z 


.i* 

2 


» = -e v ~ p 

3 * 4 


the i 


+ *, * = e 
- J 


•it- +i r -ir + „ 

-e -e + e =0 


4-36 

(d) 

(e) 

(f) 

(g) 


(h) 
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-it +i t -i* +i Y 


C j + I 77 =e 2 + e 2 + e 2 + e 2 = 2 


+1 


t -it 


e 2 + e 2 


4/ 


=4 cos — (j) 

2 


Hence 


4 U/ y 1 . 

-F=+-cos ^ +/j j+ J n(i/ - ja ) + constant 


J 


Similarly 


-dG 


dz) 


= l± 


a j 2 

< v+ u j) + 17117 




= 4 COS J +Zj. + 


J 




and 


- G = 4 cos + Zj + V ^ in(^ - constant 

2 - j J 


From relations (a) to (d), inclusive. 


'7TJJ 

v/ 


+ 1 (F + G) 


IT v 

w 


= -i(F - G) 


Thus finally 


^ C.S t-icosV^B 

ti; _ ' y ' -> 


V 


(i /-y )(^-v ) 




(k) 


(i) 


(m) 


( n ) 

(°) 


(3.9) 


ZT- = -2^ - 7)cos ^ - isinY i n -Tsin f (3 .10) 


w 


Note that according to Equations (1.44) and (1.45), -w Q cos ^ = u Q and 


-w Q sin^K = v Q where u Q and v Q denote values of u and v on the top side of 


the wing outside the Mach cone. The choice of the constants of integration 
may be justified in the following way: The relations given above for u Q 
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and v Q are in particular valid at the point whose V -coordinate is i. At 
that point Equations (3.9) and (3.10) are also valid. The first term in 
Equation (3.9) gives 0, and the real part of the second term is -w Q cos^, 
as can be seen by comparison with the corresponding term in Equation (3.3) 
But this is the correct value since cos = tan (j ) . In Equation (3.10) it 
may be noticed that the first two terms are imaginary for v = i. Since sin^ 
n/i - cos^f = s/\ - tan^ (f), the constant agrees with the value from sweep- 
back theory. 

For the case of no sweepback ( </) = 0, ^ = ' rt /2), the simplified rela- 
tions are obtained: 


7T U r~ , 1 v 

— — = -i V2 (y + jj) (3.11) 

w o 

(3.12, 


From these relations it can be seen that there is flow around the 
side edge in the following sense. On the lower surface u = 0 on the wing 
outside the Mach cone. Inside the Mach cone it increases steadily and 
reaches the value + c« at the side edge. The positive sign means outward 
flow, as is physically reasonable since there is overpressure on the lower 
surface. At the edge, v (which of course is constant on the wing) jumps 
to the value + 00 (positive sign means directed upward). On the upper sur- 
face the flow is inward, into the suction region. At the side edge u = - o o 
and increases monotonically to the value zero on the Mach cone. In the 
plane of the wing outside the wing, u has the constant value zero. 

The sidewash and downwash velocities are obtained in terms of the 
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physical coordinates r and 0 and of the sweepback <f> by taking the real 
parts of Equations (3.9) and (3.10), then making use of the trigonometric 
relations in Section III-D and relations (1.37). 

Referring to Figure 13 and formula (3.9) 


(^l)(^ 4 ) 

where 

- arg (* / - y l)(* / - y 4 )" 


TTu , 1 . 0/ . 0 

— = 2 (N - -) cos - ' sin - + cos f ' arg 


(3.13) 


may be obtained from relation (f ) in Section III-D . Hence the sidewash 
velocity inside the Mach circle is 


-TTu 


(— 3 (1 + tan^ ) , 


sin ^ “ tan 


arc cos 


1 - r (1 - cos 9 + tan <p) 

[l + r^(tan^$ - sin^G) - 2 r cos 0 tan^j 2 


(3. 13') 


In the plane of the wing the formula for the sidewash is 


t > 0 : u 


t < 0: 


-2 |— r^- ( x + tan j ■ 


l+t(2+tan^) 

tan^ arc cos ! _ t . tan ^, (3.13") 


That u = 0 off the wing in the plane of the wing follows directly from the 
discussion at the end of Section I-J. 

Referring to Figure 13 and formula (3.10) 
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? = -2 (N -jj> cos t ■ cos f + sin W.rg [<^1^1 
2 I (^-^)(^-V 4 ) 


^ ) (1 + tan </> )j cos ® - sint^(^ 2 +^ 3 )-(^ 1 +^ 4 )J + IT r (3.14) 


From relation (a) in Section III-D 


tan( 4 » 2 +«^ 3 )= 2 

(cos t + N ) + 


9 if/ 0 

2 sin ? (cos Z + N cos 2 ) 


2N cos -r • cos ~ 
^ 2 


tan(* 1+ ^> 4 )= 2 7 (cos g- - N cos 2 j 

(cost + N 2 ) - 2N cos J ' cos ^ 


tan |( 4> 2 + 4> 3 ) - (^ + ^ 4 )j = 


^ 0 

4N (1 - N^)sin 2 cos 2 

4N 2 cos 2 | - (1 + 2N 2 cos <j/ + N 4 ) 


Zrft” 1 ^ 1 - tan ^) j 


0 

cos 2 


7 0 1 

2 r cos — “ r tan 9 - 1 


Hence the downwash velocity inside the Mach circle is 


+ (1 - tam<t>) \ arc tan 


2 [r(l - r)(l - tan f)J cos j 


1 + r tan <p - Z\ 


(3.14') 


For no sweepback, 9 = 0 (and hence 9 = ^/ 2 ), Equations (3.13 1 ) and (3.14 1 ) 


reduce to 
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(b) € -Plane (b) * -Plane * 

FIGURE 14. PLANFORM Ila FIGURE 15. PLANFORM lib 


F. Planforms Ila and lib 

These planforms differ from planform Ilab by the fact that the side 
edge is yawed ( j5 t 0; Cf. Figs. 14 and 15). As pointed out in Sections 
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I-F and I-L, the solution may be obtained from that for planform Ilab by a 
Lorentz transformation. One may use either a transformation of the veloc- 
ity potential or a Lorentz transformation of the w-function. However, be- 
cause of the different boundary conditions for leading and trailing edges, 
neither of these methods will give the complete solution for both types of 
planforms; some correction functions will be needed. For this reason both 
methods will be investigated. But first some preliminary formulas will be 
derived. Consider a wing as given in Figure 14 or 15. Such a wing may be 
unyawed by the oblique transformation (Cf. formula 1.17) 

P— >P* 


P = (x lt X 2 , X,), P* = 


x^ - bx^ 


x 3 - bx l 


v 1 - b 2 ' X2, ✓n^b z ‘ 


(3.17a) 


b = tan f?> 

The corresponding homographic transformation is (Cf. Section I-L) 

1 ~V\ - b 2 


, On ^ " B 

e ~* ?i = r^Be’ B = Ts ( b ) = 


(3.17b) 


Denote the sweepback of the original yawed wing by 
transformed wing by 0°. 

Then by Equation (1.50*) 



and that of the 


tan <j)° 


tan <j)k - b 
1 - b tan<J) b 


(3.13a) 


In the complex plane Equation (3.18a) will correspond to the formula: 
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if* 


if° e* -B ,,, o , xt 

e-,=e ^Tj.re = “ b.cosV = tan <J> 

1 - Be 1 ' 


1 


(3.19b) 


Let Wq and w° be the values of w on top of the wing outside the Mach 

bob 

cone in the yawed and the unyawed cases, respectively. Define u Q) u^, v q , 
v° in a similar fashion. Then by the general formulas (1.44): 

b b jb © c>,o 

u„ = - w tan© , u = - w„ tan © 

O n r o O r 


b 

w = 
o 


-v 


o 

-Vo 


( 3 . 20 ) 


J 1 - tan^ 0° ° s/l - tan Z 0 


2 /Ao 


Using the fact that the inverse of Equation (3.18) is obtained by 
substituting -b for b (Cf. Eq. 1.20), one obtains from Equation (1.51) 


v/r 


tan Z <p° = 


'l - bVl - tan 2 0 b 

1 - b tan 


(3.21a) 


/ 2 j.b VI - b 2 v/l - tan 2 <t>° 

1 - tan <p = — — ~rZ 

' 1 + b tan p° 


(3.21b) 


If one also assumes that v° = v^ , then Equations (3.20) and (3.21) imply 


W o° - b 2 

^ b - tan</P - -2B - (1 + B^) cos W ' 


1 - B 


2 b (3.22a) 

-ltbtanf _ -(1 + B ) + 2B cost 

tantp/T^T^ ~ (1 - B 2 ) cost 6 
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w 




\/l - b 2 1 - B 2 

l + btan<£° ”fTl B^n^B^osV^° 

,b i b 

- 1 ~ b tan (/> 1 + B - 2 B cosV 

\/l - b 2 1 - B 2 


(3.22b) 


Lorentz Transformation of the Velocity Potential; Solution for 
Planform Ila 


This is the first method mentioned in Section I-L. Denote by U°,V°, 
and W^he three functions defined by Equations (3.9), (3.10), and (3.3), 
respectively. According to Equation (1.57) one may form three new func- 
tions U b , V b , and W as follows ( V denoting ^ ~ rather than 1 /?): 

“V) - 7Tp- u ° w - 


, 1 ° 
+ — 

IT o 


•a + b 2 ) cos r - zb 


1 - B 


>-')(>-'- V. 


Z)\ 3 ) 


(3.23a) 


-w y 0 . vg 

V (fc) = V°ft) = ~^jr 2 cos ~T (u- y ) + — 


(^-^3)(> / -^4) + 


(3.23b) 


wb<f) = fr^r o^rrp- w ’<?>=TiTF‘“^r 


\V 


IT 

o b 


2B cosr + 1 +B^ (*'-*' 1 )(*'-*' 4 ) 

lin (^ 2 )(^ 3 ) 


1 - B 


(3.23c) 


Since v Q = v Q when the transformation is applied to the velocity 
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potential (Cf. Eqs. 1.21c and 1.57b). Equation (3.22) is valid. With the 


aid of Equations (1.56), (3.20), and (3.22), Equations (3.23) may be sim- 
plified to 


b _ 


u 


uV) = -tT- 


V b (e) = ir 


W b (e) = If 




c 2( y " » ) + i J n 


Z n 3 J 

(^-^l)(^-^ ) 

(^-^ 3 )(^-^ 4 ) 


, . ^ 4 ) 

-)+iin (v-^)(^-^ 3 ) 


+ IT 


(3. 23a') 
(3.23b 1 ) 
(3.23c') 


Here the constants may be given in many forms; e.g., 

1 


c, = 


y 2(1 + tan £1 m - b tan0 b )(l + tan0 b ) . 


* b + tan (fP 


1 - b 


( 1 + b)tan 0 


c = / 2 - / 2(1 - b tan0 b ) 

\/ 1 - tan (tP v (1 + bj(~"l - tan0 b ) 


1 - tan ^ V (1 + b)( 1 - tan0 b ) 

bV^2(l + tan^f) _ b /l(\ + tan</) b )(l - b' tan 0 b ) 
“ 1 + b tan0° 1 + b v 1 - b 


(3.24a) 

(3.24b) 

(3.24c) 


'3 1 + b tan 0^ 

The real velocity components corresponding to Equation (3.23) may be 
easily evaluated in physical coordinates. The velocity on the original wing 
is related to that on the unyawed wing by Equation (1.21c) and the velocity 
on the latter wing is given by Equations (3.8* ), (3.13 11 ), and (3.14*). It 


should be noted that according to Equation (1.21c), v is the same at 
corresponding points. Hence assuming w^ to be the same for both 
wings, O’ is also the same. Since the sweepback angles differ, w Q also 
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takes on different values on the two wings according to Equation (3.20). 

It follows from Section I that the three functions given by Equations 
(3.23) and (3.24) represent a conical flow field. However, one has to in- 
vestigate whether they actually satisfy the boundary conditions for the 
wings in Figure 16. It follows immediately that all conditions on V are 
satisfied and that in the plane of the wing w* = 0 on the wing and w = 0 
off the wing. By construction these properties are transferred from the 
unyawed wing to the yawed wing. There remains to be checked the condition 
at the side edge. Formula (3.23) shows that for B f 0, W b has here a pole 
of order \ because of the term I/ 2 A A comparison with Equation (1.13 1 ) 
yields the result (Cf. also the discussion preceding Eq. 1.46)’: 

Formulas (3.23) and (3.24) yield the correct solution 

for planform Ha (leading side edge) but not for lib. (3.25) 

A solution for planform lib will be obtained from the second method 
of applying the oblique transformation. 

E. Lorentz Transformation of the W -function; W for Planform lib 

In accordance with Equation (1.58) one generates a function W^ b (e) 
from the function W° defined in Equation (3.3) by 


WiV) = w°ty) 

» = t ~ B 

L l - Be 


(3.26) 


As before the planform of the unyawed wing is given by Equation 
(3.18) or (3.19). Relations (3.20) and (3.21) are also valid. However, 


1 
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Equation (3.22) fails since the condition v° = v^ has to be replaced by 

w° = w*\ Thus the angle of attack and v on the wing are different in the 
00 

b o b 

two cases. If v q is given, then v q and w^ are determined by 


o b 

w = - w - 
o o 


1 - tan 2 0° \J 1 - tan 2 0 b 


(3.27) 


which follows from Equation (1.44) and w° = . 

According to Equation (3.21) one may write Equation (3.27) as 


1 + b tanC 


1 - b tan0 b n _ b 2 


(3.27') 


If these values of v Q and w q are used in determining W°, then the 
function defined by Equation (3.26) satisfies the boundary conditions on 

j d^yb , 

the Mach circle for both planforms Ha and lib. For B T 0, (— ) ~r 0. 

der e =0 

Hence by Section I-L, there is a logarithmic discontinuity in the downwash 
at the origin, and Equation (3.2 6) is good only for planform lib. Further- 
more, Equation (3,26) gives zero perturbation pressure at the edge (T[ = 0). 
This is a second reason why it gives the correct boundary condition for 
planform lib but not for planform Ila. 

In summarizing Section III-H above, 

Wjk as defined in Equation (3.26) satisfies the boundary 
conditions for planform lib. However, for planform Ila 
it fails for two reasons: First, there is not a pole in 

w at the side edge, and second, v has the correct value 
only for -1 <. t <0 and changes discontinuously at t = 0. (3.28) 
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One might, however, approach the solution of planform Ila in the 
following way. Start with the function defined in Equation (3,26). Then 
try to add a correction function which removes the difficulties mentioned 
in condition (3.2 8) without spoiling the boundary conditions which already 
are satisfied. This correction function must thus satisfy the following 
four conditions: ^ 

1. w = 0 on unit circle. 

2. w* - 0 on the wing, 

3. W is regular inside the Mach circle except for a 

simple pole in v at V - 0. ^ (3*29) 

4. v is zero for -1 < t < 0; at t = 0 it has a dis- 
continuity which is the negative of that resulting 

b 

from Wj . 

Because of condition 1, W may be extended by reflection and is hence 
regular on the whole v -plane except for the simple pole \/p at the origin 
and the reflected pole V at V - 00 . This fact, together with condition 2, 
leads to a function of the form 


W = k W 2 b = k i ( »+ J-) 


(3.30) 


where k is a real number and 


V 


e - b 


1 - Bt- 


In order to satisfy condition 4, k will have to be determined from 
the equation 
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IE k 


d(W 1 b + k W b ) = 0 (3.31) 

e = 0 

If this equation is solved for k, one obtains (Cf. Eqs. 3.23c 1 and 


3.24c) 


w 


b 


b 



c 3 


w b = Wjh + k W 2 b 


(3.32) 


where the function W is the one derived in Section III-G. 

b b 

Returning to planform lib, one sees that Wj + c W_> , c being any 
real constant, satisfies all boundary conditions except that for the side 
edge (1.13 1 ). Hence for this planform there are infinitely many solutions 
unless a Kutta condition is imposed at the trailing edge. Thus Equation 
(1.13 4 ) is actually an essential boundary condition. In planform Ila the 
situation is different. Here the value of k is determined by the downwash, 
and the existence of a pole at the edge is proved to be necessary. 

If poles of higher order are allowed, one might also consider the 
correction functions 


-n n 

W = ci( V + i/ ) ( 3 . 33 ) 

where n is an odd integer. These functions are excluded on the grounds 
that they yield infinite lift on a finite area (Cf. Section II-C). 

I. Solution W = i + (1/^)3 ; Dowrwash and Sidewash for Planform lib 
The function W ^ = i [v + (1 ) ] has already been discussed in this 

section. This function will appear in many connections and thus needs to 
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be studied separately here. In computing the corresponding sidewash and 
downwash (U^ and V^^), one has to distinguish between two cases: 

< 0 and f5 ^ 0. 

In the first, where ^ 0, several properties of the flow field may 

be determined a priori . The solution represents a flat wing of zero 
thickness and of planform lib. Although the angle of attack is everywhere 
zero, still the wing is lifting with an infinite pressure peak at the 
edge. On the Mach cone and outside the Mach cone all components of the 
perturbation velocity vanish. As pointed out in Section I-J, the sidewash 
u is = zero for 0< t <1 and = constant = u Q for b < t <0. This condi- 
tion represents a vortex sheet behind the lifting wing. Since the lift is 
positive and u Q is the value on the upper side of this sheet, u q -c 0. The 
corresponding value on the lower side is -u Q . On the upper surface of the 
wing u increases from - 00 to zero. At the edge the upwash is infinite and 
there is also a logarithmic infinity in V at t = 0 corresponding to the 
logarithmic jump of U. 


b b 

Actual computation of U 2 and bears out these predictions. One 


obtains the following trio of functions: 

U b = ELiL , y+ 1) + 


-iAn 






B 

+TT 1 

1 «l 




(3.34a) 


v 2 b - - 1 # (. 4) ♦ in ~ 7 fj* (3.34b) 

(v+Z-B)^- L_) 

W 2 b = i(^+r) v/ ' B (3.34c) 
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Here 


B < 0, V = 


/* -B 
1 - Bfc 


For evaluating Equation (3.34) in real coordinates, the following for- 
mulas may be used: 

B 2 +l 1_ lV (1-BHHB) 2 H-b (3.35) 

2B b ’ 2B ‘ b 4B /-B b ^ _2b 

t - b 

Also on the wing, putting t* = ^ ^ 


i(P+?) = i 


y it*[ 


1 


y I T*1 


1 - 1t*I J zji - lt*l ) 

/ M J |t*| 


■/ 


2(1 + t*) _ / 2(1 - bXL + t) 
-t* / b - t 


Hence 




(3-36) 


where 


-1 < t < b 


As may be seen by comparing Equations (3.23c), (3.2 8), and (3.22b), 


one may write the solution Wj for planform lib as 


W, 


W b - c W 2 b 


(3.37a) 


where 


c = 


TT 


4B 
1 - 


cos 


r 
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Hence the sidwash and downwash for planform lib are obtained from 
Equation (3.34) and the formulas 


U 


b 

1 




(3.37b) 



(3.37c) 


The second case, fj} > 0, also represents a wing of zero thickness. 
However, since b > 0, the planform is Ila and the ray t = 0 is on the 
wing. The angle of attack is zero for - 1 < t < 0 and jumps discontinu- 
ously to another constant value between 0-C t < b. Wings with such dis- 
continuities will be discussed in Section VI. 

If the sidewash and downwash are computed, one obtains the following 
trio of functions corresponding to Equations (3.34) 


U 2 b = -i Sja / i.\ . (1 - B)(1 ♦ B) Z U 

2B ' 4BS~B (u+iy/B)(u- L_) 

v/TT ’ 


b 1 - B 1 


2B 




(1 - B)(l + B; 


4B V~B~ 


-ifn 


(y+ i>/B)( 

(y- i \/b)( v 


(y+i V / B)( y + -jL ) 




(3.38a) 


(3.3 8b) 


W 2 b = i(*+^) 


(3.38c) 


As before, in Equation (3.35), 


(1 - B)(l + B) 2 
4B v/b” 



1 - b 
2b 
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Note that the second terms in the expressions for sidewash and down 
wash differ from the corresponding terms in Equation (3.34). The real 
part of the second term of Equation (3.3 8a) will now be evaluated for 
points on the wing: 

* = i J |T*| , T* = < 0 (a) 


Hence putting | T*| = R* 

(v- i\/B)(v+ —~) (n/R*" 

VB 


(i/+ i s/B)(p - — L 
B 


Here 


/5)(>/5* ^ 


) 


VR* . 1 - B 


1 - R* V / B~ 


( yiT + Sb)(Sr? -h 1 + 

1 - r^"". i-b ^ 5 1 ■ E * 

vtm b 


/' - t* . 1 

y t *+i 


^~5" 


1 - B 
/E- 

(b) 


and 


t* = 


t - b 
1 - bt 


(c) 


-t* 


b-t 


t* + 1 = (t+l)(l-b) 


(d) 


Expression (b) may also be transformed into 

.12 

\ /b - t - v b(l + t) 

yoTbjT J {e) 

Hence on the top surface of the wing (-1 < t <. b) the following for 
mulas result from Equation (3.3 8): 


Page 10Z 




(3.39c) 


J. Planform III 

Consider a flat lifting delta wing with both edges inside the Mach 
cone. To begin with, leave unspecified whether one edge is trailing 
(planform Illb) or whether both are leading (planform Ilia) and what the 
nature of the flow at the edges is. The boundary conditions are shown in 
Figure 16. 



FIGURE 16. PLANFORM III 
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The Riemann surface of the function W can easily be obtained by re- 
flection (Cf. Sections II-B and III-B). The fact that w = 0 on the unit 
circle means in geometrical language that W maps the unit circle in the 
6 -plane on the imaginary axis in the W -plane. Then by analytic continua- 
tion one can define W outside the unit circle in such a way that points 
which are images with respect to the unit circle are mapped on points 
which are images with respect to the imaginary axis in the W -plane. The 
image of e in the unit circle is 1/ e and the image of w + iw* in the 
imaginary axis is - w + iw*. Hence the rule for extension is in algebraic 
formulation (Cf. Eq. 2.9) 

W(l/?) = - W(e) (3.40) 

Before extension, W was analytic inside the unit circle except across the 
slit B 2 B 2 and except for possible poles at Bj and - By extension W will 
then be analytic in the outside of the unit circle except for the slit 
(1/Bi to + 00 , - 00 to I/B 2 ), which is the image of the slit inside the cir- 
cle. Hence W is defined on the whole Riemann sphere ( fc- -plane and point at 
infinity) except for the two slits mentioned. By Equation (1.46 ), w* = 0 
on the slit B 2 ^ both on the top side and the bottom side. By Equation 
(3.40), this condition is then also true for the reflected slit. Since 
w* = 0, one may again use analytic extension for continuing the function 
across the slits. The rule for this extension will be 

W(e) = W(7 T (3.41) 

According to this rule w = Re(W) is symmetric with respect to X-axis. But 
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from considerations in Section I it is known that in the part corresponding 
to physical reality w is antisymmetric with respect to the wing. Hence 
continuation across the slits will not give the values of W that the orig- 
inal part of the function assumed but will lead to a new sheet of the 
Riemann surface The situation is very similar to that of planform Ilab 
except that only one slit had to be considered in planform Ilab. It may 
again be concluded that the Riemann surface has only two sheets because 
w* = 0 on both sides of both slits . Hence the four possible reflections 
(in the upper and lower parts of the two slits) yield the same result, and 
two successive reflections give back to the original function (Cf . plan- 
form Ilab) . Thus the Riemann surface consists of two spheres joined along 
the two slits and — o« — - — . This condition is shown in Figure 17 

B 1 B 2 

where some closed curves have been drawn to illustrate the nature of 
the branch points . 


Y 



FIGURE 17. RIEMANN SURFACE FOR PLANFORM III 
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i£k 

This is a well-known Riemann surface (Cf. pp. 229 and 244 ff. of Ref. 18; 
also Part III, Chap. 9.1 of Ref. 19). It has the connectivity of a torus, 
and any function which is meromorphic on this surface (i.e., is analytic 
except for isolated poles) may be expressed as a rational function of the 
two quantities e and the radical 

/> = j (* -B 1 )(e-B 2 )(fr-g-)(c--^-) (3.42) 

Furthermore, for any meromorphic function on a compact Riemann sur 
face, the number of zeros equals the number of poles; when the poles and 
zeros are known, the rational function of /° and 6 can be written imme- 
diately. From the process of extension just described, it follows that 
the only possible poles are at the end points of the slits. 

First consider a basic solution which has a simple pole at B 2 and a 
zero at B^. The connection with physical condition at the two edges will 
be discussed later in this section. From Equation (3.40) one concludes 
that I/B 2 must be a simple pole and l/Bj a zero of the function W. Since 
no other poles are possible and the number of zeros already equals the 
number of poles, it may be concluded that there are no other zeros either. 
The rational function of /° and 6- which has the correct zeros and poles is 
easily seen to be 

(e- B 1 )(e-^-) 

W = k ^ 1 (3.43) 

where k is a real constant. The constant k has to be real but is otherwise 
arbitrary. Then for real values of Equation (3.43)gives the correct 
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values: w* = 0 on the slits and w = 0 between the slits. Furthermore, 
w = 0 on the unit circle. Hence Equation (3.43) satisfies all boundary- 
conditions except Equation (1.4 1 ), which will be used in determining k. 

Note that a simple pole does not have the form 


e- B 


but 




This condition exists because B£ is a branch point, and a neighborhood of 
B 2 on the Riemann surface cannot be described by the variable e; instead 
one needs a uniformizing variable (Cf. pp. 214 ff. of Ref. 18) like >/ e - B 2 . 

By a similar reasoning one sees that (e -Bj) represents a zero of order 2. 
However, in Equation (3 .43) %/ €• - Bj appears in the denominator. If this 
factor is divided out, remains in the numerator and represents a 

simple zero. To show the poles and zeros, write 


W=k 


-(e- BjMBjfe - 1) 


(3.43*) 


(e - b 2 )(b 2 e - l) 

where k is a real number. 

Next look for solutions with simple poles at Bj and B^. Such a 
solution is found by adding to Equation (3.43 1 ) a similar solution with Bj 
and B^ interchanged: 


W = k, 


'-(* - B 1 )(B 1 e - 1) 
(e - B 2 )(B 2 e - 1) 


+ k- 


'-(«■- B 2 )(B 2 e _1 ) 

(e - BjHBjfr- l) 


(3.44) 


where kj and k 2 are real numbers. Evidently Equation (3.44) satisfies the 
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boundary conditions for W since each of the terms does and the conditions 
are homogeneous: Each of the terms has the property that its real term 
vanishes on the unit circle and its imaginary term vanishes on the real 
axis between and B^. Hence the sum of the two terms has the same prop- 
erty. Furthermore, since there are poles of the right order both at B^ 
and there remains only one thing to be checked, namely, the absence of 
a logarithmic singularity at fc = 0 if this point is on the wing (planform 
Ilia). This can be achieved by choosing kj and k 2 in the correct ratio so 
that (^-^) vanishes (Cf. Section I-J and the use of the correction 
function in Section III-H). Then 



Hence for planform Ilia, Equation (3.44) must be of the form 


W = k ' 


b F5 - tEH . -b FE ~ b 2)(b 2 ^i) 

Z J ( e - B 2^ B Z 6 ~ ^ j (*- BjKBjS-- 1) 


(3.45) 


where k 1 is a real number. 


For the unyawed delta wing j3 ^ and B^ = -B^ =B. Then 

[b& 2 - (1 + B 2 ) £ + b] ± [-Bfe 2 - (1 + B Z ) fe - b] 


W = k'B 


\/(t- 2 " B^)(B 2 t 2 - 1) 


- 2k 1 B 2 (1 + e 2 ) 

- B^)(B 2 t Z -li 


(c) 
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y 


Hence for the unyawed case 


B (V 2 + 1) 


W = k 


\/(e 2 - B 2 )(B 2 6 2 - 1) 


(3.46) 


where k“ = -2k’B is a real number. 

This is the formula obtained by Stewart in Reference 16. It may also 
be obtained directly by the following reasoning. The Riemann surface for 
W is as described in connection with Figure 17. There are four poles, 
namely, at + B » and + l/fi and these are the only singularities. Thus the 
function may be written down as a rational function of e and/°(Eq. 3.42) 
as soon as the four zeros are found. Now w = 0 on unit circle. Further- 
more, for the unyawed wing the imaginary axis is a line of symmetry. Hence 
dw/3x - 0 here and thus also 3w*/3y = 0. But w* = 0 at e = 0. Hence w* = 0 
on the whole imaginary axis. Thus at Hh i, w = w* = 0 and also W = w + iw* = 0. 
But these are actually four zeros, since the surface is double-sheeted and 
the points 6- = + i occur on both sheets. Thus the information about the 
zeros and the poles is complete, and this fact yields Equation (3.46) im- 
mediately. Comparing the results with the boundary condition (1.13*) for 
the edges, one reaches the conclusion that Equation (3.43*) is the solution 
for planform Illb if B^ corresponds to the trailing edge and (Eq. 3.45) is the 
solution for planform Ilia. Both formulas contain an undetermined con- 
stant k which will have to be determined from the angle of attack (Cf. 

Eqs. 3.50 to 3.54, inclusive). 

As for planform lib, one obtains the result that if a pole is allowed 
at a trailing edge there are infinitely many solutions for planform Illb, 
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because the angle of attack gives just one relation between kj and k 2 in 
Equation (3.44). Thus this formula represents a one-parameter family of 
solutions for planform Illb, all of which have poles at the trailing edge 
except when k 2 = 0 (Cf. Ref. 16 for a discussion on higher-order 
s ingular itie s ) . 

To evaluate w in physical coordinates in the plane of the wing, ob- 
serve that 


(t-BjKBj* -1) J7 


f 2 ,l B 1 + 1 


ZB) 


1 - L. 

t bj 


(e- B 2 )(B 2 e- 1) £ - t 1 . B ?, 2 + 1 B 2 "l_ _ J_ 


F, < d > 


2 * 


Hence for the wing with one edge trailing 


w = k 


2B- 



(3.47) 


For the wing with two leading edges 

w = -k* yZ-B^lT, 



(3.48) 


Finally, for the symmetrical case (b^ = - b ^ - b), Equation (3.48) 
reduces to 


= -k'B 


J 


*-‘E> 


1 - (T f 


(3.49) 


To complete the solution the constant k has to be determined as a 
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function of the angle of attack. One way of finding k is to compute the 
downwash from the pressure and to compare it with the known values on the 
wing. Calculations have been carried out for the symmetrical case by 
Stewart (Cf. Ref. 16). For the yawed case the constants were first found 
in Reference 20. 

For planform IHb the constant k is given by (see Ref. 20) 



<* w 


l°Z 


2 E(C) + CjDf^ [c,(-Dj)] 


(3.50) 


where 


A = 


_J ~ b l b 2 + / 1 ’ b l i / T ^ 


: i =7> - c 


2 )(1 - b 2 2 ) 


D = 


V b 2 


1 + bjb 2 - b,^ 1 - b 2 2 


d 1 = y 1 - 

E(C) = the elliptic integral of the first kind 
fl 0 |c,(- D^)J = the elliptic integral of the second kind 

■ / 

(1 )J( 1 - A Z )(1 - C 2 A 2 ) 

The corresponding formula for planform Ilia is 
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k ' _ -g w oo 1 


1 - E(C) 


(3.51) 


where C is defined as in (3.50). 

In the symmetrical case Equation (3.51) reduces to 


aw • 2B 
k" = -2k' B = — 9 *— 
1 + B 2 


1 

T 7 

1 - B 2 " 

E 

1 + B* 


aw b 

^ 00 


E( v / r^b 7 ) 


(3.52; 


From Equations (1.43a), (3.49), and (3.52) it follows easily that the 
sidewash on the wing has the simple form 


-ff^U/b) 


u = 


E(7f- b z )/(l - {t/hY 


(3.53) 


With the aid of this formula and the rule for oblique transformation 
(Eq. 1.21), the constant for the yawed case of planform Ilia (Eq. 3.51) 
may be deduced from the unyawed case (Eq. 3.53). 

The downwash is determined from 


dV 

de 


; . ... . (1 + -B 2 : 

2 2 

> (i-ej 

* U w OO * 

E 

1 - B 2 
1 + B^ 


(e 2 - b 2 )(b 2 t- 2 - 1) 

3/2 


IV. WINGS OF CONSTANT-LIFT DISTRIBUTION 


(3.54) 


The wings considered in this section are lifting wings of zero thick- 
ness. In this respect they are similar to the wings discussed in Section 

♦ 

III. Thus the local angle of attack is the same at corresponding points 
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on the top surface and on the bottom surface, and v is symmetric and u and 
w are antisymmetric (Cf. Section I-B). However, in Section III a constant 
value of v was prescribed on the wing since the shape of the wing was 
given. In the present case the lift of the wing is prescribed and hence w 
is given. It will always be assumed that the lift distribution is con- 
stant. Thus on the top surface w = constant = w Q and on the bottom sur- 
face w = - w Q . For the purpose of obtaining solutions the present case is 
thus much closer to Section II than to Section III, because, in the present 
case, w is prescribed in the whole plane of the wing and also of course on 
the Mach cone. Thus there are sufficient boundary conditions for one per- 
turbation velocity. As for the symmetrical wings of given shape (Cf. Sec- 
tion II), this fact makes the solution straightforward. Another similarity 
to Section II is that interference between two wings of constant lift in 
the same plane does not affect the boundary conditions, for there is only 
downwash and sidewash interference but no pressure interference between 
two wings of the constant-lift type. 

These general characteristics simplify the solution in two ways. 
First, because of the last-mentioned property it is enough to obtain the 
solution for some very simple planforms and then obtain the other plan- 
forms by simple superposition as in Section II. By making a slight modi- 
fication of the procedure in Section II planform II will be chosen for the 
basic solutions. Planforms I and III are then easily obtained by superpo- 
sition. As a matter of fact if, e.g., solutions for planform lib have 
been found, one may obtain solutions for planform Ha by subtracting from 
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the two-dimensional case. Condition (l.B) is not of importance for the de- 
duction of the analytic formulas since for the constant-lift wings there 
is no difference in principle between a leading and a trailing side edge. 

But such a subtraction is very useful for actual numerical computations . 

The second simplification is that for solving problems in this sec- 
tion one may merely borrow the solutions from Section II. More precisely, 
if a certain function f(f-) is the correct solution for V for a symmetrical 
wing with v = + a prescribed on the wing, then the same function f(e) is a 
correct solution for W for a wing of the same planform with w = + a on the 
wing. Actually W(e) will be obtained for a wing of planform Ilab by this 
procedure; W for the general planform II is then obtained by a homographic 
transformation. This transformation is completely trivial for wings with 
prescribed lift distribution, and none of the complications of Section III 
arises. However, because of the special use of the wings in this section 
(Cf. Section IV-B), the interest is mainly in the downwash and the side- 
wash, which, of course, may not be taken immediately from Section II. 

Thus, the program is as follows: From Section II one obtains imme- 

diately the function representing W for planform II. From this function V 
and U are computed with the aid of Equation (1.42), and v and u are then 
evaluated in terms of physical coordinates. To evaluate w separately is 
of course superfluous since this function may be had immediately from 
Section II . 

A . Planform II 

This planform will be treated in detail because of its importance for 
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downwash calculations. The boundary conditions are shown in Figure 18. 

Since the lift distribution is constant on the wing, w is also constant on 

the wing by Equation (1.3a). It is assumed that the value of w on the top 

side of the wing is w . Hence it is -w on the bottom side. 

o o 



FIGURE 18. WING OF CONSTANT -LI FT DISTRIBUTION, PLANFORM II 


As in Section III-F the following abbreviations have been introduced: 


b = tan j3 , B = 

\W — 

fi l :e ' " e l 


->A 


, cos 1 ^ a tan d> 

b r 


(4,1) 


* 
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As before, move the point B, representing the edge, to the origin by 
the transformation 


*1 = 


e - b 
l - Be 


(4.2) 


The singularities £ ^ and are thereby moved to ^ and respectively, 

where 






e ,-b 


B£ 


and 


^2 = 1 - 


Be, 


(4.3) 


In the 7£-plane the boundary conditions look as shown in Figure 19. 



FIGURE 19. BOUNDARY CONDITIONS IN THE ^ -PLANE 


As is seen, there are three logarithmic singularities inside the unit cir- 


cle, namely, at ^ = 0, Figure 19 shows paths consisting of small 


half circles around the singularities. The paths are inside the circle 


and do not cross the slit on the negative real axis . If ^ varies counter- 


W 'V- 

clockwise on these paths, — increases by 1 at 7 = 0 and decreases by 1 at 

w o 
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"*7. and Y[ . A function which has these singularities (Cf. Ref. 4) is 


iw 



in + in (n - )£ 2 ) -i ni£ 


(4.4) 


On the real axis of the \ -plane the real parts of the first two terms 

represent angles of equal magnitude but opposite sign. Hence the only term 

“ iw o 0 

that gives a contribution to w = Re(W) is — — which gives the correct 

values w = +w q , -w q , or 0. 

For the applications, the components u and v rather than w are of 
interest. Explicit formulas for these components will be derived from 
Equations (4.4) and (1.42). As in Section III-E the functions F and G are 
introduced by 

_ zJL dW 

iw e- ( a ) 

o 

and 


-7T 6r 

dG = dW 
o 

The following auxiliary formulas will be needed: 


(b) 


dix _ (1 - Be-) - ( e -B)(-B) 1 - B 2 

de (1 - Bt) 2 (1 - Be) 2 


1 


e-B 

l-Bt 


6 1 - B 
1 - Bfe, 


= (1 - Be)(l - Be,) 

*1 - b 2 (* - €l ) 


(1 - Bfc)(l - Bej) 
(1 - B 2 )(e- ej 


(c) 


(d) 
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dW - iw 0 
d e T 


1 - Be j A 1 - Be 2 


(1 - B^) 


(l-Be)(e-e 1 ) ( 1 - Ber)( e - e 2 ) (e - B)(l - Be) 

2 


1 - Be 


* ' (1 - Be)(e -ej " e x 


1 - Bf, 11 1 B‘ 

_L . - — . +• 


£ + £ £ _ & 

1 1 


1 - B«e 


( e ) 


(f) 


1 - B 


1 - B 1 


£ (e- B)(l - Be ) ' Be T B(e - B) T 1 - BC- 


(g) 


<? 0 - BfeQ 


(1 - B e)(e - e-j) 


e -e T 1 - Be 


(h) 


£ (1 - B ) 


1 


{£ - B)(l - Be) " e - B T 1 - Be 


(i) 


Hence integrating formulas (a) and (b) and neglecting constants of 
integration: 


- F = 


1 ~ B€ 1 in (+ 6 - in (6 -<*) - Bin (1 - Be) 
*1 1 


1 “ Be, 1 

in e + T in (6 - f ) - B in (1 - Be) 

^ ■> t 


2 

1-B 2 


ine + 4in(e -B) -Bin (1-Be ) 


(j) 


Since 


Page 
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- G 


= e i in (e - - - in (1 - Be) + in (e -e 2 ) -gin(l - Be) 


B in (6- B) - ^in (1 - Be) 


wC = T ( F + G ) + c 
o Z 1 


(k) 


«) 
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-F + G 


+ c 2 


(m) 


where and c are constants of integration, then 


-2E?,_ f - g + C J _± 
lw o 1 V 1 


- B&2 1 - B 

L. + 


-) 


fn e 


+ <TJ +(L+ c- 2 )in(ff 2 ) 

+ B)in (1 - B*) - (g+B)J?n(tf- B) + C j 


1 + B 

= ( b Re*-j)j^n£+ 2Re(6^)-f n(( - 6 j)(^ ~ 


-(B + —)ln(6- B)(l - Be) + 


( n ) 


TTU 

w 


= 1 


Re e, - 


1 + B* 


1 2B 


Jne - iRe(^) J. n(<S -^((f - e^) 


+ I (B+ B )in( ^“ B)(1 " + C 1 


(4.5) 


_ . 1 - B«: 1 - Be 1 - B 2 

27rV =-F +G= (-— J- ! + jin* 


+ (e 2 + (ej -e 2 )in(£-e 2 ) 

+ (“ -B)in(l-Bf)+(-i+B)in(e-B)+c 
B B 2 

or 
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The most convenient way to evaluate the constants and c^ from the 
boundary conditions is first to find the real parts of V and U. Figure 20 


gives for a point e = R • e 


sin - R sin8 
an 1 ~ cos^ - R cos0 


- tan 0£ = sin + R sin0 
cos^ - R COS0 
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Hence 




— tan — 


tan + tan 0 2 
1 - tanGj tan0^ 


- 2 R sin 0 
cos ft - R cos 0 

sin^ V* - R 2 sin 2 e 

* + (cos'f' - R cos 0)^ 


-2R sinQ (cos^ - R cos 0) 

1 + R 2 (cos 2 0 - sin 2 0) - 2R cos0 cos^ 


(P) 


Similarly 

tan | Re £i X n e - e^jj = tan(0} * 0 2 ) 


2 sin ft 

cos ft - R cos 0 
R 2 sin 2 0 - sin^ ft 

i + ; 2 

(cosy - R cos 0) 


2 sin ft (-COS ft*R cos 0) 

1- R^ - 2 cos 2 V42R cos^cos© 


(q) 


To evaluate c^ and c 2 , put e = 1; i.e., R = 1, cos0 = 1. All the 
terms in Equation (4.5) have vanishing real parts, i.e., Cj = 0. All the 
terms in Equation (4.6) vanish except the term computed in formula (q) 
above. The right-hand side of formula (q) reduces to 


2 sin ft (-cost ± i) sin ft y/l - tan 2 ft 
- 2 cos^ ft +2 cost cost' tan ft 


Hence 


c 2 ~ “ (^^l) 


1 -er 


Re ( i i n j“T~ ) 




, ? J . *•« 

tan c <f> arc tan r 

r tan 


tan 2 ft 




(r) 


Introducing these values and some standard simplifications gives 
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JZk 


= i ^ tan ^ - ctn^Jin ( - i tan fi J n(c -6 j)(£ ~£ 2 ) 
+ i ctnj6in(€ - B)(l - Be) 


(4.7a) 


and 



(ctnjj - tan 4> )ine + i 




+ ( 


B Z -1 

2B 


) in 


6- - B 
1 - Bfc 



arc tan 


Vl - tan^/ 
tan <f> 


(4.7b) 


To complete the evaluation of the real parts the third term in 
Equation (4.7a) should also be considered. With reference to Figure 21 

Re (iJ n(e - B)(l - Be)] = - (6 1 + 0 2 ) + n TT 


where 


n = + 1 for B 0 
= 0 for B < 0 


Y 



FIGURE 21. f -PLANE 
The angles and 0^ are determined by 
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tan 0 i = 


R sin 0 


* R cos© - B 


„ R sin 0 

tan 0 -> = r 

2 R cos 0 - L 
B 


tan ( 0 J + 0 2 ) = 


R sin 0 fz R cos 0 -(B + B )J 

R^(cos 0 - sin 0) - R cos 0 (B H — ) + 1 

B 


(s) 


Hence 

ifu / , 1 v ^ . - 2R sin 0 (tan 4> - R cos 0) 

— = (-tan <f> + r~) 0 + tan <f> arc tan 7 5 5 L 

0 D 1- R^sin 0-cos^0)-ZR cos0tan <p 


-to 


arc tan 


1 _ 

2 R sin 0 (R cos 0 - b ) 

— 7 2 2 “ 7 ? 

R^(cos 0 - sin 0) - 2“ cos0 

b 


+ nir) 

+ 1 7 


(4.8a) 


where n - 0 for b < 0, n = -1 for b > 0. Actually n TT = - argument of 

(-B). The values of arc tan and n have to be chosen so that u = 0 for 

* = 1 . 

jj £• - B _ , x 

To find the real part of X n ^ ^ in Equation (4.7b) it may be ob- 

served that if 


then 


and 


| e - B | r Rj and 1 - B * = R. 


Rj 2 = R 2 + B 2 - 2BR cos 0 


R 2 2 = 1 + B 2 R 2 - 2 BR cos 0 


Hence 
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Re Unyf^)=inR 1 - in R z = £ i n 


2 2 

H + B -2BR cos0 


2BR cos 0 


The expression for v then becomes 


= (ctnjS - tan£) X nR+^\/l-tan ^ 


2 v/l-tan 2 cft (R cos0 - tan</ >) 
arc tan j_p 2_2 tan^ +2R tan <f> cos 0 


- arc tan 




tan 2 (ft 


tan <ft 


B 2 -! 


in 


R 2 4 B 2 - 2BR cos 0 
1 + B 2 R 2 “ 2BR cos 0 


4B , n^t>2 _ 7 q id a (4.8b) 

Formulas (4.8a) and (4.8b) will now be specialized to points in the 
plane of the wing. For such points 0 is either 0 , + TT (top side), or - 7T 
(bottom side). This value has to be substituted in Equations (4.8a) and 
(4.8b), or Equation (4.7) may be used. Care has to be exercised in choos- 
ing the values of the arc tan function. A reference to Figures 20 and 21 
is here recommended. 

The following calculations apply to the top side of the plane of the 
wing. Because of antisymmetry, u takes on the corresponding values of op- 
posite sign on the bottom side. 

The expression for u becomes 

tan </> + ctn^J argt - ctn ^ jarg (e - B) + arg (1 - B*)J 


vT o = (-tan <f> + ctn ( 2 >) nj - n 2 ctn £3 

where 

nj = 1 for t < 0 . ri£ - 1 for t < b 

= O for t > 0. =0 for t y b 

The corresponding formula for v reads 


(4.9a) 
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= (ctnfj - tan^>).£ n |t|+ 1-tan 2 ^ 


arc tan 


,r. 


tan 4> (T - tan <ft ) 


1 - T 2 -2 tan 2 ^+2T tanf) 


arc tan 




1 - tan ^ 
tan <p 


+B 2 - 1 , 

B - T 

+ 2B in 

BT-1 

„ 



(4.9b) 


Planform II will now be specialized in four different ways by assuming 
simple values of the angles <{> and yS : (1) <^> = 0, (2) j3 = - ^/4, (3) <f) = 0, 
fl > = -^4, and (4) j3 = 0. In the first case (<p = 0), Equations (4.7) and 

(4.8) reduce to 


^*-£[ in *- Jln < € " B) < 1 ~ Be >] 


(4.10a) 


u A 

=-ctnp 


TTV _ B+l j t-i. B-l i f*B T 

*o 2B -* n e -1 ^ n 6+i + 2B 1-Bfe 2 

2R sin 0 (R cos 0 - ctn ft ) 


-Q + arc tan^ 


2 (cos 2 0-sin^0)-R cos 0 ctn 


1+nTT 

+ 1 / 


(4.10b) 


(4.11a) 


TTV 1 fl n i j. 2 R COS 0 

rn- = — in R + arc tan 

o b 1 - R 2 

b 2 -i 2 n R 2 ± b 2 ~ 2BR cos 8 _ X 

4B 1 + B 2 R 2 - 2BR cos0 2 
For points in the plane of the wing, there is a further reduction: 


(4.11b) 


u 

W- 


= (n^ - n^) ctn p 


* 


n^ = 1 for t <. 0, iij = 0 for t > 0 
n^ = 1 for t<^b, n£ = 0 for t >b 


= t- i n j T | + arc tan 


2T 
1 - T 2 


B Z -1 

2B 



TT 

2 


(4.12a) 


(4.12b) 
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In the second case 


(j3= - ’V 4 ) 


one obtains 


ITU 

^ = i (tan +l)ine - i tan <pJl n(e - 6^(6 -(f^J-Zi in(l+*) (4.13a) 


= - (1 + tan^)in£ + i \/l - tan ^ in^—r^ 


/, ^2^ x v 1 - tan^ 

“V 1 - tan 9 arc tan I 

tan 9 

-2R sin 0 (tan^ -R cos 0 ) 

- (tand> +l) 0 +tan <b ’arc tan j j- . ^ 7 " 

v r ' r 1+R (cos 0-sm^0)-2R cos 0 tan p 


+ 2 arc tan 


R sin 0 
1 + R cos 0 


tt„ . I 2Vl-tan 1 2 ^ (R cos0-tan <f> ) 

= -(l+tan<i>) J! nR + yl - tan 6 arc tan . . 2 jl * 1 Q 

w 0 v T/ ' l-R"-2 tan^ +2R tan <p cos 0 


- arc tan- tan ^ j 
For points in the plane of the wing: 


'1 - tan" 


= - (tan <f) + l)n^ 


= 1 for t < 0 , nj = 0 for t > 0 


= “(1 + tan^>) Jin T 


2 V 1 - tan^ <b ( T - tan 


Wl - tan f arc tan x . x 2 . 2 tan 2 ^ + 2T tan ^ 


arc tan 


1 - tan 


u 1 tun j j 1 

tan <f> 

In the third case (<£ = 0, /3 = - 7 T /4) one obtains 
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= iinc- - 2 iJL n(l + f) 
o 


(4.16a) 

ss. 1/.::; -i 

o 


(4.16b) 

It U Q n , ^ S ^ n ® 


(4.17a) 

W” = ~ 0 + 2 arc tan l t) Q 

w o 1 + K cos 0 


■tTv g , 2R cos 0 

w- = --<nK + arc tan — •% — 

O 1 - R* 

r 

2 

(4.17b) 


For points in the plane of the wing: 

— - - 1 for t < 0, u = 0 for t >0 
w o 

TTv n 1 „| 2T T T 

= -Jin T/+ arc tan z - — 

o 1 I 1 _ 2 

In the fourth special case ( jS - 0) the following limits have to be 
evaluated: 


(4.18a) 

(4.18b) 


lim IE ilZ BH 1 ~ B*Mne _ lim ln(fr- B)(l - B^)-Jn^ 


B-*0 


B — >0 


2B 


1 r(5B in<t ' B ><* - B,t >] B=0 


1 / 

7( e + 7) 


N 


lim 

B-*0 


1 „ e - B . lfd n & " B 

~2B ( Jn 1 - Bf J ^ JLn 1 - Be 


1,1 

ji-z *«) (b) 


B=0 


Then Equations (4.7a) and (4.7b) reduce to 


^ = i tanf> ine - i tan* In(t-f 1 )(t*6 2 ) - jf) 


(4.19a) 
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= - tan <f> 1 


£. ^ ^ \^1 ■ +fl«^ P n 4- ~ I • 


6 

tan^ i n T^e“ + i( F * 6 ) 


•/ 


-y 1-tan ft arc tan 


v/TT 


tan ft 


(4.19b) 


tan ft 


The corresponding real parts are most easily found by comparing Equa- 
tions (4.19a) and (4.19b) with Equations '4.7a), (4.7b), (4.8a), and (4.8b). 

- 2R sin9(tanft - R cos8) 

"wj = - 9 tantf> + tan<^ arc tan R 2 ( s i n 2 e - C os 2 0) - 2R tan <p cos 0 


112 ® [ /o . 

+ 7 l R tJ 


(4.20a) 


TTv 


= - tan ft i nR 


♦/7 


tan^ft 


arc tan 


\/T 


tan^ft 


tan 


* 


v/i- 


arc tan 


tan ft (R cosQ - tan ft> ) 


1 - R -2tan^ft> + 2R tanftcos© 


c °s0 .1 

+ - 3 — (o - R) 


(4.20b) 


B. Note on Applications 

Probably the most important application of the wings of constant-lift 
distribution is in finding the induced flow field behind a wing of finite 
chord. This application was indicated by Busemann (Cf. Ref. 4) and is 
discussed in detail in Reference 2. 

V. MIXED TYPE OF LIFTING WING 
The problems studied in this section will be those in which the lift 
distribution is prescribed over one region of the wing and in which the 
shape is prescribed for the rest of the wing. Only a special case will be 
treated with the following properties: In one region the lift distribution 
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is constant, w = + w Q , and the rest of the wing is a flat plate at zero 
angle of attack. The importance of such solutions is that the interfer- 
ence of a lifting element is determined (Cf. Section V-D). In Section IV 
wings of constant-lift distribution were discussed. In the plane of the 
wing such wings induce a downwash field but no pressure field off the wing. 
However, if a flat plate at zero angle of attack is introduced into this 
downwash field, an interference effect results. Since the downwash has to 
be zero at the flat plate, an additional flow field is generated which 
cancels the downwash due to the lifting element at the plate, and at 
the same time yields lift on the flat plate. Thus, though actually the 
flat plate is at zero angle of attack, it has a variable effective angle 
of attack because of the downwash field. The resulting composite flow 
field will be determined in this section, and the lift induced by a flat 
lifting element will be discussed in Section VI. 

Solutions in this section will be straightforward, because of the 
previously developed technique, in particular the use of the flatness con- 
dition (1.46 ) in combination with analytic extension by reflection and 
also the use of the Lorentz transformation. 

The previous classification of wings by planforms is no longer suffi- 
cient. Only certain representative examples of the various possible cases 
will be treated, 

A. Example I: Pressure Distribution 

The most important case is probably the one shown in Figure 22. This 
wing consists of two regions. Region I is flat and at zero angle of attack 
(v =0), and region II has a constant-lift distribution w = w q . The more 
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general case <7+0 on region I can be obtained easily by simply adding a 
solution from Section III which has the same planform as region I and which 
has a constant angle of attack +- 0. Such a wing would have w = 0 on region 
II and hence its superposition would not spoil the boundary condition: 
w = t w Q on region II. It is also evident that when the flat plate (region 
I) is at zero angle of attack, the sweepback angle of the leading edge of 
region I is irrelevant as long as this edge is outside the Mach cone; in 
this case any perturbation velocity existing in region I is induced by re- 
gion II and there are no disturbances outside the Mach cone. 

The boundary conditions in the e -plane are shown in Figure 23. 

The fact that Ct = 0 on region I and on the wing outside the Mach cone 

has as a consequence that w = 0 on the Mach circle. Furthermore, since v 
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FIGURE 23. BOUNDARY CONDITIONS IN THE e -PLANE, EXAMPLE I 
is constant on line OP (Cf. Fig. 23), it is assumed as before that w* = 0 
here by Equation (1.46*). As in Section III the function W is extended by 
reflection in OP. Assume, to begin with, that there is a function W which 
is analytic inside the unit circle except for the slits OP and OB. By re- 
flection in OP, W is extended to a two-sheeted surface consisting of two 
unit circles. The point 0 is then a branch point, and passage from one 
sheet to the other is across the line OP. On the lower sheet a slit OB is 
then obtained similar to the one on the upper surface. As before, the 
surface is made single sheeted by the transformation 

v = s/7 ( 5 . 1 ) 

The boundary conditions in the v -plane are shown in Figure 24. 

This is a simple type of boundary value problem. The narrow unyawed 
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FIGURE 24. BOUNDARY CONDITIONS IN THE -PLANE, EXAMPLE I 
delta wing (planform III) has. these boundary values for w in the €• -plane 
if it has constant lift (Cf. Section IV) and for v in the 6 -plane if it is 
nonlifting and has a symmetrical wedge profile (Cf. Section II). There are 
two logarithmic singularities at + fa. Since w = 0 on Mach circle, two 
additional singularities are obtained at +1//B by reflection. By the 
standard procedure previously used in this report, the strength of the 
logarithmic singularities can be determined, and when the singularities 
are known, the function may be written down immediately 


(^- /b)(i - /i^) 


w -l w Q 0 

W = TT 1 n (*»+ V r B)(l + >/! v) 


(5.2) 


This discussion of the singularities does not determine whether a 
constant is to be added to the right-hand side of Equation (5.2). That 
Equation (5.2) actually is the correct expression if suitable values of 
the multivalued logarithm function are taken will be seen from the evalu- 
ation of w in terms of variables in physical space (Cf. proof of Eq. 5.3). 
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For the applications, one is primarily interested in the lift distri- 
bution in region I. Since it is at zero angle of attack, it would have 
zero lift by itself. However, the presence of the lifting surface, region 
II, induces a lift in region I. This region corresponds to the negative 
real axis in the e- -plane and to the imaginary axis in the a/ -plane. A ray 
in region I is characterized by the angle X , which is negative. 


t = tan X. 



(a) 


and also 

n = v^r 


(b) 


The ray X then corresponds to the point Ni (or -Ni) in the 1 / -plane. 



* FIGURE 25. V -PLANE: EXAMPLE I 

From Equation (5.2) it follows that the value of w at the point Ni is 

» 
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given by 


V - 4*1 + 4>z - 4$ -44 ' ir ~ 0 1 + e 2 ■ w 


(c) 


where the notation is as explained in Figure 25. Equation (c) is still 
true if the point P is on the unit circle, in which case it gives w = 0, 
the correct boundary condition. It also gives the correct value w = w Q on 
interval (- /b, V^B) which corresponds to region II on the wing. 

From Equation (c) one may evaluate w as a function of X on region I. 


tan 


N * 2 "/B N 


tan 


A -u Q + 1 

e i + e 2 n + n7b 


N(B + 1) 


1- n2 . 

, N' 


Vb (n 2 -i) 


COS (0^ 4 ©2) 


12 -n 2 


1 + 


(N^-l)^ B 
N 2 (B+1)2 


(N 2 -1) 2 B 


(d) 

( e ) 
(f) 


But by Equation (1.3 7c) 

iSiilf-. t§i±lj . 2 +2 = Z( L + 1) 2( bH 

B \ 2B/ b ’ v b 7 

and 

(N Z -1) 2 (- T - l) 2 /Tf+l\ _ zlzl 

N 2 - T \ -2T/ * 


(g) 


00 


Introducing Equations (g) and (h) into Equation (f): 
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cos(0i + 0^) 


1+-L_b±i 

1 + t b 


_t b + 1 

1 + t b 


b + t + 2bt 


b - t 


(i) 


Using the formula (arc cos x) - rr = arc cos (-x). Equations (c) and (i) 
result in 


w 1 b + t + 2bt 

-w-^-arc cos— p-g— 

t = tan t (5.3) 

For the top side of region I, arc cos is to be taken between 0 and + -rr. 
The boundary between regions I and II corresponds to t = 0. Here w = 

( w o/tT ) arc cos (b/-b) = w 0 , which checks. Where the Mach cone cuts region 
I, t = -I. Then w/w q = arc cos (b - 1 - 2b) / (- 1 - b) = 0. This value 
also checks since region II cannot induce any lift outside the Mach cone. 

B. Example I; Sidewash and Downwash 

U and V can be calculated most readily following the procedure in 
Section III-E. Introducing V = '/i- in Equations (1.42) 


Let 


dU _ 
d^ 

i 1 2 1 dW 

" + i7r> 5T 

<«) 

dV _ 
d^ 

1 , 2 _ * dW 

2 v 2 

(b) 

dF _ 

_ _]£_ 1 dW_ 

( c ) 

d v 

iw i/ ^ d is 
o 


i) 

U| * 
^ Iti 

-J-.S dw 

iw 0 d 

(d) 


Differentiating Equation (5.2) 
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IZk 


dW 
d u 


iw. 


TT 


- 4- 


* '-Jk ^ " + ^s- 


Expanding -^ 5 - and w^ dW into partial fractions 
» c dv 


(e) 


dF /B + l" | _1 

di; = " 2 ( VB J » 




B(iJ- \/b) _L B(v+ /B) iJ + L_ 

Vb Vb^ 


1 / - /b 


Integrating Equations (f) and (g) 


f'-i) y+vff 4 + i) 


(f) 


(g) 


F = 2^1 V + 


/ B+l \ 1 

W5/ 


+ constant 


G = Z 1 — — — 1^4- 


Wb ) 


+ constant 


-In (is- v/B) + Bin(l - uJ B)- - in(i/+\ZB)-Bln(l+Xv€j 
B B 

(h) 

Bln(y-\/B) + “In(l -V + v€j- -/n(l+^yB) 

B B 


(i) 


Then the complex sidewash and downwash are 

jr — = y(F+ G)+ constant 


w Q 2 


7 f Y— = |(-F+G) + constant 


Hence 

-rr 

7 T 


il = 7b±iYL 4if^4 1 

W ° Vn/bA 7 Lb/ (V+ v6)(i vS) 

r jg^K. 1). 

w o V V® A v / \ 2 B / (v»+ /B)(l -J/ VB) 


(j) 

00 


constant (5.4) 
(5.5) 

+ constant 


Page 136 



NACA TN I685 

JPl 


Progress Report No. 4-36 


The evaluation of the real sidewash and downwash (Cf. Eqs. 5.6 and 
5.7) shows that both constants of integration in Equations ( 5 . 4 ) and (5.5) 
must be zero in order to satisfy the boundary conditions. 


Of primary interest is the sidewash in the plane of the wing. Since 



The constant Cj = 0 in order to satisfy the condition U = 0attf=-1 (or 

R = 1, 0 = + tt ); 


Transforming Equation (5.4 1 ) to physical coordinates gives, with the 
aid of Equations (1.36), (1.37c), and (1. 3 7 i). 


Tu 

VV 


2(1 + b) 
b 


^(1 - 1 1 1 ) 

\*\ 


sm j 



(I) 


For the sidewash in the plane of the wing three different cases have 
to be considered, namely, (- 1 ^ t < 0 ), (0 < t < b), and (b < t < 1 ). 

In Region I, sin0/2 is = 1, - 1 .< t < 0, and hence j t | = - t. The 
last term in Equation (i) is given by Equation (5.3). Hence the sidewash 
is 


Fu / (1 + b)(l + t) 1 / b + t + 2 bt ^ ( 5 

*r 0 = Z V (Hbt) * b arccos l ' t - b ) ’ 

for 

(-1 < t < 0 ). 

In region II, sin0/2 = 0, the first term of Equation (5.4 ) vanishes 
and Re(TTW/w 0 ) = TT . Hence the sidewash is constant. 
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Tf u _ 7 T (5.6b) 

w o " " b 


for 


(0 < t <b) 

Finally for b < t < 1 it follows from Section 1-J that u = 0. 

Taking the real part of Equation (5.5), the downwash is found to be 


TTv ( B + lV l 

W °V v^/l 



V - + u /b) 


(v+ V®)(1 -v /b) 


+ C- 


(m) 


The constant c^ = 0 in order to satisfy the condition v = 0 at e = -1 


(R = 1, 0 = + 7T , V = +i). 

The downwash can now be obtained in the plane of the wing: In region 
I, 0 = + TT and > is imaginary. Thus v = 0, which checks the boundary con- 
ditions. In region II, G = 0 and y = Then 


(v - /B)(l + WS) 


( /T - v/B)(l + /bt) 


(jy + vS)(l -v v/b) 


( n/T + /B)(l - v/BT) 




t 2 +i b 2 +i t - i : 

2T + 2B + Vr~ 

l " B -2 
7 b 



t 2 +i b 2 +i 


I ZT " 2B I 

Then with the aid of the formulas for the Tschaplygin transformation 


(J - /B)(l +v n/b) 

b + t - 2 v4(l - t)b(l- b) - 2bt 


b - t 


Hence the downwash in region II is 


TTv 

VT~ 

O 


/(I + b)(l - t) n /T- b 2 


bt 


in 


b + t - 2bt - 2\/bt(l-t)(l-b) 


b - t 


(5.7) 
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The solution will of course be obtained from that for Example I by a 
Lorentz transformation. The situation is very similar to that in Section 
III (planform II). The transformation of the velocity potential may be 
used. The functions U, V, and W for Example II are then obtained imme- 
diately from Equations (1.57), (5.2), (5.4), and (5.5) and the values of 
u, v, and w in the plane of the wing from Equations (1.21c), (5.3), (5.6), 

and (5.7). However, the constant w 0 has to be adjusted. If w 0 ° is the 

c 

constant value of w on top of region II in Example I and w Q the corre- 
sponding value for Example II, the formulas (1.21c) give the following 
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relation between these two values: 

c 1 


w 


O 

n/ 1-c2 

C = tan y, 


“ c 



b>(l-c 2 

b-c 


w o ° (5.8) 


This equation shows that w Q C actually is a constant and also that w Q 

b-c 


bv/l-c 2 


w 


in the formulas (5.2), (5.4), and (5.5) has to be taken equal to 
in c.der that the method of Lorentz transformation of the 
velocity potential give the prescribed value w Q C in region II of the yawed 
wing . 

The details of this transformation may be easily carried through as 
in Section III. However , only the other possibility will be worked out here , 
namely, the oblique transformation of the pressure function. 

The solution will be written W = Wj + W 2 where W x is given by Equa- 
tion (5.2) after unyawing the wing by a homographic transformation 

-iw Q (W-V f B*')(l-V n/b») 


W l = 


JL n 


+ v r B*)(i+ n Vis*) 


(5.9a) 


where 


-■A -B* =?f# B .C=Ts(c) 

and the correction function W£ is 0 for Y < 0 and is of the same form as 
Equation (3.30) for Y ? 0 


W ? = 


ikw O . 1 


v + r) 


(5.9b) 


where k is a real number and must be evaluated to satisfy the condition 


Page 140 



NACA TN 1685 

JPL Progress Report No. 4-36 


d(Wj + 


dfr 


= 0 


(5.10) 


|« = 0 

That this is the right correction function is shown as in Section III. 
Note that Re(W 2 ) = 0 for V real and positive. Hence the addition of W ^ 
does not violate the boundary condition that w = w Q = constant in region 
II. The evaluation of k follows: 


ilT dW 


w 


o d v 


- - 2 \/B* 


*= 0 


1 


1 


C + B* 1 + CB* 


2(1 + B) s/i (B - C)(l - CB) 
B(1 - C 2 ) 


irr dWd 


W" 


o dy 


* - k <> + V 


l€= 0 


Then 


d(W 1 + w 2 ) 


d v 


and 


iw. 




TT 


2(1 + B) y ( B - C)(l - CB) k(l + C ) 


B(1 - C 2 ) 


= 0 


k = 


~ 2C(1 + B) >/( B - C)(l - CB) 


B(1 - C )(1 + C) 

With the aid of Equation (1.3 7) this formula is easily transformed into 


(5.11) 


k = - 


c /2(1 + b)(b - c) 

b(l + c) V i - c 


(5.11' ) 


Finally the pressure function is 
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•fW- TT „ (v -VB*~)(1 -V \/B*~) 

— = — (w +w ) = - iin =zr* 7=r- 

O w o 12 ( + \/B*) (1 + ^ \/b*) 


(5.12) 


2 i C(1+BMB-C)(1-CB) . 

B(l-C 2 )(l + C) 1 v’ 


where 


* = ^ 

Vi - c fc 

The evaluation of w = in the plane of the wing in terms of 

the physical coordinates follows: 

1 1 1 — A* 

=Re ik(^i) = k — VR* sin — for 0* = 0,+ TT (5.13) 

w o V s/b* 2 " 

= +k -rf=f “/T*i for 0* = + TT 

— v / I T >5c I • “ 


= 0 for 0 * = 0 

In region I in the plane of the wing with the aid of Equation (1.3 7i) 


-= = k 


\/-T* 


(5.14) 


for -1 4. t* < 0 


-1 <. t c 


where 


1 - ct 


The expression for w^ is obtained from Equation (5.3), replacing b by b* 


and t by t*. Hence, the pressure in region I is given by the formula 
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ILw = arc cos b* + t* + 2b*t* / 2(1 + t*) 

w o t* - b* + J (-t*) 


(5.15) 


for - 1 < t* < 0 
- 1 ± t c 
where 

t* 


t - c 


b*-. 


b - c 


1 - ct ’ “l - cb 
and k = 0 for /< 0 and k is given by Equation (5.11) or (5.11*) for / > 0. 

Inserting the values for t*, b*. and k, one may rewrite Equation 
(5.15) as 


» 


TT 


W 


— = arc cos ik. ~ -Ell! + *) t (1 filil 1 b) £ c < 0 
(t - b)(l + c) 


-rr 


w 


= arc cos C M 1 + *) + (*“ 0( 1 4 b ) 
o (t - b)(l + c) 


_2c /(b - c)(l + b)(l ± t) 


b(l + c) 


c - t 


for c > 0 


(5.15'a) 


(5.15'b) 


These formulas are valid on the flat part of the wing, i.e., - 1 < t < c. 

D . Note on Applications 

The use of the mixed type of wing for obtaining solutions is illus- 
trated in Figure 27. 

The solution for the flat lifting delta wing AjBA^ is given in Sec- 
tion III-J. If one wishes to have the solution for the corresponding wing 
with the tips Aj and cut off, one has to take the basic solution 

for A 1 BA 2 and then cancel the lift in the cut-off region by infinitely 

many mixed -type wings. A typical wing of this kind would have (a) constant 
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FIGURE 27. MODIFIED NARROW DELTA WING 
lift on region D 2 E 2 F 2 (region II in Fig. 22) and (b) zero downwash to the 
left of this region (region I in Fig. 22) . Condition (b) is necessary be- 
cause the basic solution for already has the correct downwash on the 

wing D^C^BC^D^, anc * *he superposition of additional solutions should not 
lead to a violation of the boundary condition that the wing is flat. 

Another application is to wings of very low aspect ratio. The sim- 
plest example is the rectangular wing. If the aspect ratio is greater than 
unity, one obtains the solution by superposition of the solutions given in 
Section II-E and a two dimensional flow. However, if the aspect ratio is 
smaller than unity, one also needs solutions of the type studied in Sec- 
tion V-A (Cf. pp. 42 and 43 of Ref. 2). 

Many similar applications are possible (Cf. also Section VI-E) . 
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VI. WINGS WITH DISCONTINUOUS SLOPE 

This section discusses conical wings with prescribed shape which are 
piecewise flat; that is, the surface of the wing consists of conical re- 
gions where (X = constant and the value of the constant is different in 
adjacent regions. Thus CC is discontinuous on the boundary between such 
regions. 

The symmetrical case is solved by trivial superposition of the solu- 
tions given in Section II. This is the same superposition principle that 
was used in Section II. As was explained there, it is justified because 
there is no downwash interference between symmetrical wings in the same 
plane. 

The lifting case is also trivial when the upper and lower surfaces 
are independent, i.e., when all free edges are supersonic. In this case 
it is reduced to the symmetrical case as explained in Section III-A. This 
method is illustrated in Examples I and II of this section. 

The case where at least one edge is subsonic is less easy. Examples 
III, IV, and V belong to this type. Although these examples may not be 
obtained from any symmetrical case, it will be seen that some results pre- 
viously obtained in Section III will furnish part of the solution. 

It will always be assumed that there is no gap at the place of dis- 
continuity of angle of attack. By gap is meant a break in the solid sur- 
face through which the region above and below the wing are connected. In 
such a gap the pressure would have to be zero in the plane of the wing. 

If the gap is infinitely narrow, however, there is only a removable 
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singularity; i.e., the pressure would have to be zero along the ray repre- 
senting the gap, but this value would not influence the rest of the solu- 
tion. It must be remembered, though, that a gap may have finite extension 
in the vertical direction and still its projection on the plane of the 
wing have zero width. This condition is easily visualized for a deflected 
control surface (Cf. Section VI-E). However, su*fch a finite gap cannot be 
treated by the present theory which applies all boundary conditions on the 
wing to the corresponding projections in the plane of the wing. 

Although there is a great variety of wings falling under the heading 
of this section, only a few examples are treated. They are in general of 
direct or indirect practical importance and are also chosen so as to il- 
lustrate the methods. 



FIGURE 28. PLAN FORM FOR EXAMPLE I 


A. Example I 

This wing is illustrated in Figure 28. It may be reduced to a simpler 
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case. First of all, it is enough to consider the case of a symmetrical 
m wedge profile. Then angle of attack means half angle of the wedge. From 

this case one obtains the solution for the flat plate profile by the meth- 
od discussed at the beginning of Section II. Furthermore, it is sufficient 
to afssume that the angle of attack = (dj -# 2 ) = ^ ^e ^ ne 

of discontinuity and = 0 to the right of this line. Then the general case 
may be obtained by superimposing a wing of the same planform and constant 
slope, -#2 (Cf. Section M-B). 

Hence the problem has been reduced to the wing shown in Figure 29, 
but here the flat plate is immaterial. Thus there is simply the wing of 

* 



FIGURE 29. EXAMPLE I, REDUCED 
planform II whose solution was given in Section II -A, for such a wing does 
not cause any downwash in the plane of the wing outside of the wing, as 
has been pointed out repeatedly. Hence a flat plate may always be inserted 


* 
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there without changing flow conditions. 

Thus the pressure distribution on the wing in Figure 29 is given by- 
formula (2.8c) for /3 = 0, and formula (2.8c) combined with formula (2.13a) 
for 0. 

The principle used for treating Example I is evidently quite general. 
One may superimpose several nonlifting wings to form a nonlifting wing 
with prescribed slope distribution. If necessary, flat plates are then 
inserted to prevent any communication between the upper and lower surfaces. 
After they are inserted, one modifies the wing so that the lower surface 
is intact but the profile changes to a flat plate profile. Then one has a 
lifting wing whose angle of attack at any point is equal to half the wedge 
angle of the unmodified wing. The flow below the wing is the same in both 
cases, and the conditions above the wing are obtained by the general sum- 
metry principles of Section I. One more example of this principle is il- 
lustrated in Figure 30. 





Profile: Symmetrical Wedge or Flat Plate 
FIGURE 30. PLANFORM FOR EXAMPLE II 
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B . Example II 

By the general principles stated in Section VI- A, one obtains the 
solution for Example II immediately from the solution for planform III of 
Section II and Example I of the present section. 

C . Example III 

The example illustrated in Figure 31 may not be solved by the above 



Profile: Flat Plate t p is Counted Positive 
FIGURE 31. PLANFORM FOR EXAMPLE III 
methods: The case of the nonlifting wing with wedge planform is of course 
easily solved by superposition of solutions from Section HI. But since 
the upper and lower surfaces are not independent, the solution for the 
flat profile may not be deduced from the solution for the wedge 
profile . 

Fortunately, the boundary conditions for Example III are very sim- 
ilar to the conditions (3.29) for the correction function in Section III-H. 

It can easily be seen from the discussion in that section that the solution 
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for Example III is given by Equation (3.30) 


w = i^-( if + y) 

TT v >> 


( 6 . 1 ) 


where k is a real number and 


, /=■ ,, 6 - B „ 1 - v4 - b 2 , _ . ^ 

^ -V 7 ! ’ ^ t . B ~ ■ B = , b - tan fi 


To make this formula complete, one has to establish the numerical 
relation between the constant k and the discontinuity in angle of attack 
at a = 0. This relation is obtained from Equation (1.42c). 


1_ . (1 - B«) - ( * - B)(-B) 


dW _ ik , J_v . 1_ . 

J ' - ' " 2 /Jf (1 - B*) 2 


d6 TT 


(a) 


(-) 
'de ' 


6=0 



1 _ 1_ _ 1 - B 2 

2 * /^B * 1 


_k. . 1 ~ B • 1 + B 2 + 2B 
+ /b 2B 




k /z ( 1 - b) 

“ i TT V b 


•(k + D 


(b) 


In evaluating Equation (b), Equations (1.37c) and (1.3 7i) were used, 
and the undetermined sign was introduced because of the square root. The 
correct sign of k may be determined in the final formula by the rule that 
w is positive for the positive value of (X . Hence by Equation (b) and the 
discussion following Equation (1.42c) 


a 
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* 


Discontinuity in v = ft 

= +iri (- i) - il - B X‘ ± B ) 2 

2 T 2B Vb (c) 

+ TT i(- i) * fM Z b ) b + 1 

2 TT V b b 


Hence by the above rule for the determination of the sign of k 

4B \/b OC w„ 


k = + 


b / 2b 
1 + b / 1 - b 


ftVL = 


LOQ 


00 (1 - B)(l + B) 2 


( 6 . 2 ) 


if 

"=+ ym on top of wing 

In proving Equation (6.2) one should remember that 
^ 1 -i 

vnrr 

is the dominant term in Equation (6.1). 

Formulas (3.38) of Section III —I give the corresponding sidewash and 
downwash. Thus the complete solution for Example III is the following 
trio of functions 


* 


JLy = - 2 V~b q + b 2 ) 

(1 - B)(l + B) 2 


Tl V _ 2 V'eT 1 _ 
1 + B ^ V 


-v) + 


TT_W = 4B /B 

(1 - B)(l + B) 2 




(v + i-/B) (V- 

/B 

-i/„ (‘'-^>("--75) _ ^ 

(v , +i'/B)(^ + 

- 

i(^ + i) 

)/ 


Here 


(6.3a) 


(6.3b) 


(6.3c) 
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V = 


/«- B 


2 vB 1 + B‘ 


2b 


1 - Be 1 - B (1 +B)‘ J 1 - b 1 + b 


2Vb" 


2b 


4B/B~ 


2b 


1 + B y 1 + b ’ (i_b)( 1+B) 2 7!- b 1 + b 

It is easily seen that v has the correct discontinuity at 6 =0. 

The evaluation of formulas (6.3) in physical coordinates in the plane 
of the wing can be made easily from Section III— I . In particular , the 
pressure on the wing is determined by 


w = 


gw «, 2 b yr rrrr 

r l + b J b - t 


■l <t <b 


(6.4) 


D . Example IV 

This case is a generalization of the previous case as is seen from 
Figure 32. It is of course reduced to the previous case by the transformation 

^ -c 


1 - Ce 


-= 0, C 


= L ~ 7 I C — , c = tan/ (d) 



FIGURE 32 . PLANFORM FOR EXAMPLE IV 
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The solution is obtained by combining formulas (1.57) and (6.3). One has 
to substitute 0 for e and B* = (B - C) / (1 - BC) for B in Equation (6.3) 
and then form linear combinations of these expressions according to Equa- 
tion (1.57). The combination of these substitutions leaves the definition 
of V in Equation (6.3) formally invariant, for V should now be (0-B*)/(l-B*0) 
which is simply (e - B)/(l - Be), as may be seen by direct computation. A 
more intuitive way of proving this invariance is given in the following 

paragraph: 

Z 

l) is the result of two successive homographic transformations of the 
special type discussed in Section I-L,. The first one takes B into B*, and 
the second one B* into zero. Hence the composite transformation takes B 
into zero and has the form given above independent of C. 

The procedure for obtaining the solution is straightforward and re- 
sults will be given only for W : 


TTW , 1, 2C « (i J - i /B*)( 1 / + — 4 = ) 

=k,i (* + -)+ -in' n \T B *' 

^ w oo !- c (t/tiyFKy- i \ 


~ 7b* 


- B 
1 - B 


B - C 
1 - BC 


4B* n/b* C(1 -t- B* 2 ) t 1 j- C 2 

1 ~ (1 - B*)(l + B*) 2 B*(l - C 2 ) 1 - C 2 
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From Equation (6.5) one obtains the following formula for the pressure 
on the upper surface of the wing: 



One application is to find the lift of a twisted wing. The simplest 
example is a wing with rectangular planform and zero thickness where the 
slope is constant chordwise but varies spanwise. Such a wing may be ob- 
tained by superposition of wings of the type shown in Figure 29 with 
01 = 02 = /3 = 0 • However, it must be assumed that the slope of the wing 
is constant within the Mach cones from the leading edge tips. If this is 
not the case, the superposition described would give lift outside the wing. 
However, this lift may be removed by additional superposition of wings de- 
scribed in Section V. Another way of expressing this fact is to say that 
the twisted rectangular wing may be obtained by superposition of wings of 
the type shown in Figure 33. 

If the solution for the wing in Figure 29 is applied to the wing in 
Figure 33, the region ABC will be lifting. This lift may be cancelled by 
superposition of wings described in Figure 22, These wings will then change 
the lift distribution in the region ABD of Figure 33. 

The problem of the lift of a flat wing in a variable downwash field 
is of course similar to that of a twisted wing in a uniform flow field. 
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FLAT PLATE 


Finally consider the problem of finding the forces induced by the 
deflection of control surfaces. Two examples are shown in Figure 34. 



I ' E F 

FIGURE 34. VARIOUS TYPES OF CONTROL SURFACES 


It is clear that the effect of the control surface ABCD in the region 
DAjA^ is given by the solution for Example I in this section. However, 
normally this region contributes just a small part of the lift. On the 
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main part of the control surface, two-dimensional theory applies. The 

effect of the control surface EGF in the region E^FG is given by the solu- • 

tion for Example IV. In this case one is no longer dealing with a small 

correction. The main effect of deflecting the control surface is not 

given by any formula in the previous sections, but only by formula (6.6). 

It has been assumed here that a gap which might exist along AD or EG 
in Figure 34 does not influence the pressure distribution appreciably. 

Further discussions of the applications may be found in Reference 3 
where all solutions given in this section have been shown to have practi- 
cal applications to control surfaces. 
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TABLE I 



NOMENCLATURE 

Defined 



on page 

P 

pressure 

4 

P <*=» 

free-stream pressure 

4 

m 

1 

~l — 2 

4 


VM 4 - 1 


M 

free-stream Mach number 

3 


/ 2 2 


r 

V X 1 + x 2 

27 

R 

m / \ l-\/l-r^ 

Ts(r) _ r 

29 

t 

x l/ x 3 

31 

Ts(a) 

1 -v/l - a 2 

Tschaplygin transformation 

a 

29 

u 

perturbation velocity in Xj-direction (sidewash) 

4 

u* 

harmonic conjugate of u 

33 

U 

u + iu* 

35 

V 

perturbation velocity in -direction (upwash) 

4 

v* 

harmonic conjugate of v 

33 

V 

v + iv* 

33 

w 

perturbation velocity in X 3 -direction 

4 

w* 

harmonic conjugate of w 

33 

W QO 

free-stream velocity 

4 

W 

w + iw* 

33 
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x 

x 1,x 2 ,x 3 


y 

a 

e 

4 

6 
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coordinate in fc -plane = Re(e) 

coordinates in physical space 

xj in lateral direction 

x 2 in vertical direction 

X 3 in free-stream direction 

coordinate in £ -plane = lm(dr) 

angle of attack 

R • e l ® = x + iy 
i 8 

r * e = xj + lx^ 

B 

1 - B* 

X y 

arc tan ~ 

X I 


Defined 
on page 

30 


30 

6 

30 

30 
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TABLE I (Cont’d) 

Defined 
on page 


half-angle of wedge 57 

Mach angle = arc tan m 54 

79 

frefc-stream density 4 

arc tan t 31 

sweep-back angle 41 

velocity potential 4 

arc cos (tan <j>) 42 


* 


> 
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